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I. INTRODUCTION 


A. BACKGROUND 

In September 1987, due to the pressing environmental issues of ozone depletion and 
global warming, an international conference was held in Montreal, Canada. This 
conference was viewed by many as a major step toward correcting the critical 
environmental problems facing the world today. Specifically, one of the problems 
addressed was chemical products damaging the protective bubble we live in. This bubble 
or thin layer is known as the ozone layer, and is what makes the earth's atmosphere 
habitable and hence unique from any other known planet. Without this protective layer, the 
earth's climate would not support life as we know it today. Protecting this ozone layer 
from chemical destruction 1s paramount. 

A group of chemicals that deplete the ozone layer is chlorofluorocarbons (CFCs). 
CFCs are manmade chemicals of chlorine, fluorine and carbon and were developed almost 
60 years ago by General Motors. They are unique in that they have a combination of 
desirable properties not often found in chemistry; low in toxicity , non-flammable, non- 
corrosive, non-explosive, extremely stable and compatible with many other materials. 
However, it 1s this extreme stability which causes problems to the ozone layer, because 
they only break down in the upper atmosphere (stratosphere) when subjected to intense 
ultraviolet radiation. This breakdown produces chlorine which has been linked to the 
depletion of the earth's protective ozone layer. A chlorine atom destroys an ozone molecule 
(03) by reacting with it and converting it to an ordinary oxygen molecule (02) and carbon 
monoxide. More disturbing is the fact that the chlorine atom survives the ozone-destroying 


reaction and is free to repeat the process - one CFC molecule being able to destroy 100,000 


ozone molecules. Unfortunately, a typical chlorine molecule can survive for up to 100 
years, so a major problem does exist. 

Armed with these facts, 24 nations (representing the United Nations Environment 
Program (UNEP)) met and signed the Montreal Protocol in 1987. They discussed 
substances that deplete the ozone layer [Ref. 1] and called for a near-term freeze on the 
production and consumption of these substances. To be more specific, it required 
production of these chemicals to be cut back to 1986 levels followed by a two-phased 
reduction culminating in cutbacks of 50% by mid-1998; this came into effect on July 1, 
1989. In 1990, a progress meeting was held in London where UNEP delegates agreed to 
completely phase out all CFCs by the year 2000 [Ref. 2]. This was a needed, but 
worrisome change for the U.S. Navy as they use a number of different CFCs for various 
refrigeration and solvent/cleaning needs. One such need is the approximately 1850 
shipboard air conditioning (AC) plants presently using CFC-12 (in reciprocating 
compressor plants) and R-114 (in centrifugal compressor plants). These AC plants not 
only provide personal comfort, but cool vital electronic components,weapons systems and 
other auxiliary equipment. 

To comply with the Montreal Protocol and U.S. legislation, the Mechanical Systems 
Branch/(Code 2722) at David Taylor Research Center (DTRC) is pursuing research aimed 
at eliminating shipboard use of CFCs. As mentioned by Krinsky and Noel [Ref. 3] this 
research is to be completed in three phases: 

1. In the short term, to identify suitable alternative ozone-safe chemicals to directly 
replace CFC-12 and CFC-114. In order to do this, the heat transfer characteristics must be 
Similar to the existing refrigerants in place and hence the need for a database exists for 
current refrigerants (i.e. R-114, R-12, R-113) so that they can be compared to the new 


proposed refrigerants. 


2. In the longer term, to research, develop and test substitute chemicals and 
alternative technologies to replace existing CFC uses. 

3. To implement new cooling system technologies into the Fleet which do not 
depend on CFCs or their replacements. One such technology under investigation at the 
Naval Postgraduate School is thermoacoustic refrigeration, which uses sound waves to 
compress and expand harmless gas molecules (such as helium) inside a tube. 

The purpose of this thesis is to supplement DTRC's research on alternatives to CFCs 
by establishing baseline nucleate pool boiling data of R-113 from a small bundle of 
enhanced tubes representing a section of a flooded evaporator . Emphasis is placed on 
surface history affecting activation/deactivation of nucleation sites, hysteresis, and the 
affect of varying evaporator pool height. The data obtained are expected to serve as 
reference data for comparison with future proposed refrigerants from bundles of both 
smooth and enhanced tube surfaces. It is felt that what is learned from this thesis about the 
characteristics of R-113 can be applied to any general refrigerant and can help the Navy 


design future AC systems using alternative chemicals. 


B. OBJECTIVES 
The objectives of this thesis are as follows: 

1. Understand in greater detail both the convection and nucleate pool boiling 
phenomena and hysteresis effects within a small enhanced tube bundle. 

2. Vary the evaporator pool height above the tube bundle to see how this effects 
nucleation site activation and general evaporator performance. 

3. Investigate the effect of surface history on nucleation site activation/deactivation. 

4. Compare data with similar studies (i.e. same fluid or tube) done at the Naval 


Postgraduate School. 


II. LITERATURE SURVEY 


A. GENERAL INTRODUCTION 

Due to the desirability of operating many engineering devices in the nucleate boiling 
region, many types of enhanced tubes are being studied. Nishikawa et al. [Ref. 4] noted 
that there are two primary ways to increase the heat transfer performance of a boiling tube. 
The first was to decrease the wettability of the fluid in which you are immersed. They 
suggested using a Teflon coating on the boiling surface. However, this is not very 
effective for refrigerants as they are very highly wetting on any surface and fill the cavities 
of the Teflon surface. The second way proposed was to manufacture a tube with many re- 
entrant cavities, thereby trapping the vapor and keeping the nucleation sites active. These 
re-entrant cavities are designed such that if one site is activated, neighboring sites will soon 
follow and cause the entire tube to nucleate almost instantaneously due to the 
interconnecting passageways. There are now many such tubes with re-entrant cavities on 
the market. 

The two main types of reboilers analyzed today are kettle reboilers and full bundle 
boilers. The difference is described in detail in Payvar [Ref. 5]. To summarize, in the 
kettle type reboiler the tube bundle only occupies the middle lower portion of the evaporator 
leaving large areas on the sides of the bundle for fluid circulation. The full bundle boilers 
occupy the whole cross section of the lower half of the evaporator shell. As can be 
reasoned from these definitions, circulation patterns are very different for the two types of 
boilers. 

The following is a brief summary of the work done on the bundle apparatus at the 


Naval Postgraduate School: 


bundle or how surface history affects Turbo-B nucleation site activation/deactivation. In 
addition, no information was found on varying the pool height above the bundle to see how 
this affects nucleation site activation and general evaporator performance. These areas 
could be very important to the U.S. Navy as well as to industry and could cause the startup 
procedure of refrigeration units to be altered, thus emphasizing the need for the objectives 


of this thesis. 


and in understanding the natural convection and nucleate boiling regions for a single tube 


apparatus and several correlations have been adapted. However, the tube bundle opens a 


T.J. Murphy, September 1987 [Ref. 6] 

¢ Designed System 

¢ Conducted Preliminary Operation 

C.L. Anderson, June 1989 [Ref. 7] 

¢ Tested Smooth Tubes, R-113 

¢ Tested Smooth Tubes, R-114/Oil 

¢ Tested Finned Tubes (19 fpi), R-114/Oil 
N. Akcasayar, December 1989 [Ref. 8] 

¢ Tested Finned Tubes (19 fpi), R-114/Oil 
¢ Tested High Flux Tubes, R-114/Oil 

H. Eraydin, December 1990 [Ref. 9] 

¢ Tested High Flux Tubes, R-114/Oil 

¢ Tested Turbo-B Tubes, R-114/Oil 


EXPERIMENTAL AND THEORETICAL STUDIES 


In the literature survey, no papers were found on the Turbo-B tube operating in a 


In the past decade, there have been significant improvements made to tube surfaces 


new arena of topics. 


Yilmaz and Westwater [Ref. 10] studied the effect of velocity past a smooth tube and 
how it effects the heat transfer performance. They used a 6.5 mm diameter tube in R-113. 
They found that with an increase in velocity, the heat transfer performance of the tube was 
better. 

Cornwell and Schuller [Ref.11] stated that the increases in heat transfer coefficient 
found by Yilmaz and Westwater could not entirely be explained by the increase in velocity 
past the tube. Therefore, they decided to conduct a photographic study of boiling R-113 in 
a bundle at one atmosphere. Their study used 241 electrically heated tubes, 19.1 mm in 
diameter and arranged in a square in-line pitch of 25.4 mm, giving a pitch-to-diameter ratio 
of 1.33. They found that the bubbles leaving the lower tubes in the bundle impacted and 
caused a sliding motion around the upper tubes. Cornwell and Schuller attribute the 
increase in the upper tubes heat transfer performance on this impact and sliding motion of 
bubbles from lower tubes. Cornwell and Schuller also showed that the heat transfer 
coefficient was the highest for the top tube in the bundle and decreased as one moved down 
the bundle; the difference was as much as four times. Cornwell [Ref. 12] later found that 
in the nucleate boiling region, sliding bubbles and liquid forced convection could account 
for all the heat transfer in the top of the bundle. 

Chan and Shoukri [Ref. 13] studied the boiling characteristics of a small enhanced 
multitube bundle in a pool of R-113. The rectangular bundle had a vertical pitch of 23.8 
mm and a horizontal pitch of 31.75 mm. From there studies, they concluded that the 
convective contribution is highest for tubes near the top of the bundle due to the increased 
quality and fluid velocity across these tubes. They also noted at lower heat fluxes, the 
convective effects dominate whilst for higher heat fluxes, the primary mode of heat transfer 
was nucleate boiling. At high heat fluxes there are numerous active nucleation sites and a 


high bubble generation rate which decreased the influence of nsing bubbles generated from 


tubes below. They also found that the convective effect from a boiling tube below an upper 
tube in a bundle has less of an effect as it gets farther from the tube of interest. 

A study by Bergles et al. [Ref. 14] showed that if a tube is on the verge of nucleation, 
it burst into nucleation if tapped. This test was repeated many times with the same results. 
Several other experimenters using different tubes, refrigerants, and bundle arrangements, 
have also shown the influence of tube position on heat transfer performance in the bundle. 
Wallner [Ref. 15] used a smooth tube bundle with both square and triangular geometric 
arrangements giving pitch-to-diameter ratios of 1.33 to 1.5 respectively. He found a 50 
percent increase in the heat transfer coefficient over a single tube for tubes near the top of 
the bundle. 

Fujita et al. [Ref. 16] used both smooth and enhanced tubes to study the effects of 
tube position, distribution of heat flux, and system pressure changes (0.1 to 1 MPa) within 
the bundle. They used a tnangular tube arrangement with R-113. In these tests, they 
reported similar results as above with higher tubes giving improved performance. The 
enhancement in heat transfer again was believed to be caused by bubbles generated from 
the tubes below. 

Webb et al. [Ref. 17] pointed out that in air conditioning applications, the refrigerant 
typically enters the bundle with a 15 percent quality. With this entering quality into the 
bundle, it simulates the effect of lower tubes activated in the bundle and thus will improve 
the performance of the entire bundle. Furthermore, the boiling mechanism for a flooded 
evaporator is very different from the kettle type evaporator and therefore it is very hard to 
formulate and evaluate a theoretical model for a bundle. 

Leong and Cornwell [Ref. 18] looked at the same representative slice of a flooded 
evaporator that Cornwell and Schuller [Ref. 11] and Cornwell [Ref. 12] studied. They 


concluded that the heat transfer coefficients were the highest at the top of the bundle and 
decreased in performance as you went down the bundle. 

Based on the above-mentioned studies, all experimenters seem in agreement, 
independent of the fluid or tube used in the bundle, that the presence of boiling lower tubes 
Causes a convective flow in the evaporator that gives mise to an increase in the heat transfer 


coefficient on the upper tubes. 


ITI. EXPERIMENTAL APPARATUS 


A. TEST APPARATUS OVERVIEW 

The whole experimental apparatus including the auxiliary equipment and the 
evaporator/condenser test apparatus used during this investigation is shown in Figure 1. 
The following is only a general description of the whole experimental apparatus. A more 
detailed look at the condenser and evaporator is provided in Section D. Further information 
about the apparatus is provided by Murphy [Ref. 6] and Anderson [Ref. 7]. 

Looking at Figure 1, the apparatus uses three closed loops. The first loop consists of 
an 8 ton refrigeration unit located outside the laboratory which is used to cool an ethylene 
glycol/water mixture. The second loop is the ethylene glycol/water mixture flowing 
through the condenser. This flow rate through the condenser was delivered by two pumps 
which could be operated independently or together; this coolant mixture condenses the 
refrigerant vapor in the condenser and maintains system pressure. Pump number one 
provided coolant flow through the four test tubes as well as to one of the auxiliary coils. 
Pump number two provided coolant through a manifold which distributed the coolant to the 
remaining four coils in the condenser. The third loop is the evaporator/condenser itself. 
The vapor generated in the evaporator flowed upward and condensed in the condenser 


which then returned the condensate to the evaporator via gravity. 


B. AUXILIARY EQUIPMENT 
1. 8 Ton Refrigeration Unit 
This unit was used to cool the ethylene glycol/water mixture (coolant) to the 
desired temperature needed to condense the refrigerant vapor. For the present pool boiling 


experiments, the temperature control was set to maintain the sump at 5°C. This 


refrigeration unit had a cooling capacity of 28 kW (8 tons) with a 30 gpm pump providing 
circulation of the coolant to achieve the set temperature. To prevent damage of the unit, 
several safety features were provided: low temperature cutout, high pressure cutout, low 
flow cutout, and a hot gas bypass which allowed for continuous operation. For operation 
procedures, see Appendix E. 
2. Ethylene Glycol/Water Mixture 
The coolant used was a 54% by weight ethylene glycol/water mixture. This. 
was used to control the system pressure by circulating it through the condenser. 
3. Pumps 
There were two pumps available to circulate the coolant from the sump through 
the condenser. Pump number one fed four test tubes and one of the secondary condenser 
coils. Pump number two was used to feed the other four secondary condenser coils. 
Pump number one provided enough cooling for tests up to five evaporator tubes. 
However, during bundle and simulation evaporator operation, both pumps were needed at 
heat fluxes above approximately 60,000 W/m~. 
4. Flowmeters 
Four calibrated float-type flowmeters, connected to pump number one, were 
used to measure the flowrate passing through the four test condenser tubes. One additional 
flowmeter which was connected to pump number two was used to measure the total flow to 
the secondary condenser coils. In addition, each of the four secondary condenser coils had 
a globe valve to further regulate (or shut off) flow as desired. In the test tubes, the flow 


was regulated by a ball valve. 


C. DATA ACQUISITION SYSTEM/INSTRUMENTATION 
As described by Akcasayar [Ref. 8], a Hewlett Packard HP-3497A Data Acquisition 


System, HP-9125 computer and HP-701 printer were used for data acquisition, data 
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reduction and data printing respectively. Although an HP-9826 computer and HP-7470A 
plotter can be used for final graph printing, a Macintosh Classic computer using 
CricketGraph 1.3.2 was utilized. HP Basic 3.01 was used for data reduction. At the 
beginning of every set of runs, the HP-9125 computer had to be initialized by using three 
system discs if the system was turned off. As described by Anderson [Ref. 7], type-T 
copper-constantan thermocouple measurements (mvolts) were made on the HP-3497A with 
the relay multiplexer assembly equipped with thermocouple compensation. A 20 channel 
relay multiplexer card was used to measure the voltage output from voltage and amperage 
sensors. Woltage measurements were taken from separate sensors that measured the 
voltage going to the tube bundle, simulation and auxiliary heaters. The total amperage 
going through the auxiliary and simulation heaters was measured using an American 
Aerospace Control (AAC) current sensor. The currents of each instrumented tube heater 
were measured using five identical current sensors. The voltage supplied to the other active 
tubes was also measured but the current of each active tube was not (the total current for a 
pair of active tubes was measured). This was felt to be sufficient since these tubes each had 
the same power output (1000 W) as the test tube heaters and there was no apparent reason 
to monitor each active tube heat flux individually. 

Computer channel assignments for data acquisition and array assignments are given 


in Table 1. 


D. EVAPORATOR/CONDENSER 

An overall view of the evaporator/condenser is shown in Figure 2. The evaporator 
was designed to simulate a small portion of a refrigerant flooded evaporator. Front and 
side views of the evaporator are shown in Figures 3 and 4. It was fabricated from stainless 
steel plate and formed into a short cylinder 610 mm in diameter and 241 mm long. 


Electrically heated tubes were cantilever-mounted from the back wall of the evaporator to 


1] 


permit viewing along the axis of the tubes through the lower of two glass/Plexiglass 
windows mounted in the front. The glass was used on the refrigerant side in order to 
prevent crazing or cracking of the Plexiglass and the Plexiglass gave the glass added 
strength. Figure 5 is a schematic sectional view of the evaporator that shows the four kinds 
of heated tubes used in this study. The four kinds of heater tubes were as follows: 
Instrumented, Active, Auxiliary, and Simulation. Table 2 gives the power rating for the 
heaters and the number used in the evaporator. 

The electric power can be applied separately to each set of heaters using a STACO 
240V, 23.5 KVA rheostat controller shown in Figure 6. Also, the desired number of 
instrumented tubes, active tubes, auxiliary, or simulation heaters can be turned on by using 
individual circuit breakers. The four auxiliary heaters, each with a maximum rating of 4 
kW, were installed two on each side of the test bundle to maintain the liquid pool at 
saturated conditions and to provide system pressure control. The five simulation heaters, 
also capable of 4 kW each, were mounted below the test bundle in order to simulate 
additional tube rows in a larger bundle and to provide inlet vapor quality into the bottom of 
the test bundle as suggested by Webb [Ref. 17]. The auxiliary and simulation heaters had 
to be used with caution to avoid breakdown of the refrigerant. 

The tube bundle consists of instrumented, active, and dummy tubes. The location of 
each tube 1s represented by the respective letter 1, A, and D shown in Figure 5. The test 
bundle consisted of two types of heated tubes: active tubes (marked "A") which contained 
1 kW cartridge heaters, and instrumented tubes (marked"I") which, in addition to the 1 kW 
cartridge heaters, contained six wall thermocouples. In measuring boiling heat-transfer 
coefficients, great care must be exercised with the cartridge heater and temperature 
measuring instrumentation to ensure good accuracy. For example, the type of heater used 


could have an affect on heat transfer. Bergles [Ref. 19] stated that the method of heating is 
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‘immaterial except for the case of a.c. resistance heating which may cause bubble generation 
in phase with the power supply. However, more recently, Jung and Bergles [Ref. 20] 
who carried out extensive pool boiling data with R-113, concluded that the heat transfer 
coefficient of a single tube in pool boiling is not sensitive to variations in the cartridge 
heater heat flux provided that enough thermocouples are used to measure an average wall 
temperature. Prior to this investigation, the instrumented test tubes were fabricated in a 
similar way to those used by Hahne and Muller [Ref. 21] and Wanniarachchi et al. [Ref. 
22]. The exact procedure can be found in Eraydin [Ref. 9]. Figure 7 is a cross-sectional 
sketch of an instrumented tube used during this thesis, showing the construction details and 
the location of the wall thermocouples. The thermocouples were imbeded in the wall at 
different circumferential and longitudinal positions along the heated section of the tube. 
The instrumented tubes were located along the centerline of the tube bundle, forming a 
vertical in-line column. All the instrumented and active tubes were Turbo-B tubes made by 
Wolverine. These tubes were 14.15 mm in diameter and were arranged in an equilateral 
triangular pitch (i.e. centerline-to-centerline spacing) of 19.1 mm, giving a pitch-to-_ 
diameter ratio of 1.35. The bundle also contained a number of unheated dummy smooth 
tubes (marked "D") that were used to guide the two-phase mixture through the bundle. The 
dummy tubes were made from commercially available 15.9 mm in diameter smooth copper 
tubing. Two vertical baffle plates were made of aluminum to restrict circulation into and 
out of the bundle from the sides . A small open space approximately 5 mm in height was 
left between the bundle and the dummy tube rack. This space allowed some liquid to enter 
the bundle and replace the vapor being generated in the bundle. However, there was also a 
space below the dummy rack that provided the majority of the circulation. Thus, 
liquid/vapor circulation was vertically upward over the five instrumented test tubes with no 


net horizontal component. The dummy rack (Figure 8) used 12 tubes made of solid 
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aluminum rods each 15.9 mm in diameter and spaced 19 mm from centerline to centerline. 
This rack had a triangular pitch arrangement with vapor retaining plates on the side and was 
used for two purposes. This dummy rack was designed to collect all rising two-phase flow 
generated by the simulation heaters and direct it to the test section. The other purpose was 
to simulate the vapor passing though a larger bundle before reaching the instrumented 
tubes. The liquid-vapor mixture after passing through the bundle was separated when it 
reached the pool surface. The vapor flowed toward the condenser and the liquid was 
circulated back toward the bottom of the bundle. 

The condenser included four test tubes (each 2 m in length) in a vertical in-line 
column and five auxiliary copper coils used to regulate the pressure in the system. It was 
designed to permit independent condensation studies of small in-line tube bundles, using 
the evaporator as a source of vapor. All these tubes and coils were cooled by a refrigerated 
mixture of water and ethylene glycol. Each of the test tubes and auxiliary coils could be 
operated separately using ball or globe valves respectively to regulate the desired flow rate. 
The vapor that was produced in the evaporator was guided up toward the top of the 
condenser by a vapor shroud as shown in Figure 9. The vapor was then directed down 
over the condenser tubes and coils where it condensed. Once condensed, the force of 
gravity returned the condensate back to the evaporator. Observation windows in both the 


condenser and the evaporator allowed for easy viewing of condensation and boiling. 


E. GEOMETRY OF ENHANCED SURFACE USED 
1. Turbo-B Tube 
The Turbo-B tube, manufactured by Wolverine Tube Inc., has a enhanced 
surface geometry. The exterior boiling enhancement is made by raising integral low fins, 
Cutting diagonally across these fins, and then rolling the fins to compress them to form 


mushroom-like pedestals [Ref. 23]. This process forms numerous re-entraint 
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passageways. Figure 10 shows the surface of the tube at 25 times its actual size. The tube 
is currently available in copper, cupro-nickel, and low carbon steel. The relative 
dimensions of the tube used in the present experiments are: 

Tube material - Copper 

Outside diameter = 14.15 mm 

Inside diameter = 12.7 mm 


Enhanced surface length = 203.2 mm 
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Figure 8. Photograph of Dummy Tube Rack 
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Figure 10. Close-up View of Turbo-B Tube Surface (25° Times) 
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TABLE 1. COMPUTER/DATA ACQUISITION ASSIGNMENT 


Thermocouple Description 
Vapor 1-Top of Condenser | 00 i T(0) 

























[Liquid 2-Top of bundle | _04 | __T1a) 
| Liquid 3-Bottom of bundle { 05] (5) 
Tube No? OCS) 
[Tube 1,No.6 +445 SC«dYCSSCiA)—C*S 
ube 2.Not dS) 
[Tube 2,No.3 Sd BCC) 
Tube 2,No.4P Cd) 
[Tube 2,No.5 Ss SC*dSCOSC A) 
[Tube 2,No.6 +41 Sid SS) 
Tube 3.No.t Sid Sidi) 
ube 3,No2 Sid SCC) 
[Tube 3,No.4 PSCC) 
[Tube 3,No.5 SiS SC*dCSCSC*i PX) 
ube 4.NotSCdSC—iBSCdYSCSCi 2) 
Tube 4,No.2 iY SSCS) 


Tube 4,No. 3 | GO T(26) 
1(27) 
Tube 4,No. 5 (28) 











| __ Tube 5,No.4_ | 4 CT(80) 
Tube 5,No.2 0 | 65 CT(31) 
| Tube 5No.3 | G82) 


Tube 5,No. 4 T(33) 
Tube 5,No. 5 cc T(34) 
Tube 5,No. 6 [eS Ge | T(35) 
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TABLE 1. COMPUTER/DATA ACQUISITION ASSIGNMENT (CONT.) 










Amp(7) 
Amp(8) 
Amp(9) 



















Auxiliary Heaters Amp(10) 
Simulation Heaters Amp(11) 








| Simulation Heaters | 28 | Voit) 






TABLE 2. EVAPORATOR HEATERS 


Heater Type Number Power Rating per Heater 


|Instrumented Tube Heaters | 5 | 1000W_ 
[| Auxiliary Heaters [4 T ODOW 
| Simulation Heaters | 5 | OW 
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IV. EXPERIMENTAL PROCEDURES 


A. REMOVAL OF TUBE BUNDLE FROM EVAPORATOR 

Before removing the tube bundle from the evaporator, it was ensured that the 
evaporator was free of any refrigerant and was at atmospheric pressure. The tube bundle 
support block bolts were then loosened and removed. The tube bundle was then removed 
from the back of the boiler. Care was taken when removing the tube bundle as to not 
damage the thermocouple or heater wires. It was easier to disconnect the thermocouple 
wires from the data acquisition board and the heater wires from the power panel instead of 
leaving them connected to the tubes when it was necessary to work on or clean the tubes in 
the bundle. The tube bundle consisted of five instrumented heated Turbo-B tubes, 10 
active heated Turbo-B tubes, 18 dummy smooth tubes and two dummy Turbo-B tubes. 
These two dummy Turbo-B tubes were equipped with heaters and could be connected if 


desired to study various other bundle configurations (see Anderson [Ref. 7]). 


B. BUNDLE DISASSEMBLY 

When the bundle needed to be disassembled, it was ensured that there was a clean 
working space and a container to put small parts in during the procedure. The first process 
was to remove the perspex plate attached to the aluminum baffle plate by four screws. This 
plate was modified such that the eight dummy smooth tubes surrounding the instrumented 
tubes extended through the plate (with ptfe plugs which held the tubes in the proper pitch- 
to-diameter alignment). This alteration became necessary since at high heat fluxes, the 
smooth tubes (which were screwed into the tube support block) loosened from the wall and 
rattled around. The ten screws on each aluminum baffle plate were then removed. The 


plates were then pulled off the bundle. The two outside smooth dummy tubes remained on 
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the tube bundle support block as they were countersunk into the block. The three inner 
smooth tubes could be removed from the aluminum block as they were attached only by the 
screws already removed. The 12 smooth tubes were then unscrewed from the tube bundle 
support block as seen in Figure 11. Care was taken in noting the order in which all the 
tubes were removed as it saved time in reassembly. The bundle was engraved to help in 
this process. With these tubes and baffle plates removed, only the instrumented and active 
heater tubes remained. These tubes were removed by loosening the outer O-ring 
compression plate; this allowed the tubes to be pulled from the tube support block. Care 
was taken when pulling the tube from the tube support block so as to not break the 


thermocouple or heater wires. 


C. SYSTEM CLEAN-UP 

If the system was contaminated it had to be thoroughly cleaned. In order to do this, 
the entire apparatus had to be taken completely apart and cleaned in the following manner. 

After removal of the R-113 from the evaporator and system at atmospheric pressure, 
all electrical connections to the bundle were disconnected and the front-viewing glass 
windows were removed. The tube bundle was then be removed as described above. The 
dummy tube rack was then taken out; at all times, the system was ventilated with a fan. 

Having removed the tubes from the tube bundle, they were individually washed with 
warm water, rinsed and then wiped down with acetone. The smooth tubes were cleaned 
with Copper Brite to remove any tarnish. They were wiped down with warm water and 
acetone. The same procedure was followed for the Turbo-B tubes except they were not 
cleaned with Copper Brite. During the cleaning process, a soft bristled toothbrush was 
used to ensure the enhanced surface was cleaned properly, exercising care not to interfere 
with the tube surface. The storage tank, if thought to be contaminated also was drained, 


vented, opened and cleaned using the same procedure. 
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D. SYSTEM REFIT 
1. Installation of tube bundle in evaporator 

Prior to installing the evaporator tubes and support block, the bundle and 
evaporator was thoroughly wiped down with acetone. Once the evaporator tubes were 
installed in the support block, and prior to tightening the tube bundle support block nuts, 
the bundle was guided into the evaporator section. After the bundle was in position, it was 
then checked that the dummy rack was still properly positioned below the bundle and that 
the vapor thermocouple positions were still 1.75 cm above the bundle. The support block 
nuts were then tightened equally on opposite sides to give equal compression of the gasket. 
If the front window was removed during maintenance, the window was put back on very 
carefully with small equal torques applied to each nut, in turn, in a clockwise direction 
around the outer ring support. After the window was in place, each tube (which extended 
through the outer O-ring compression plate) was tapped forward so as to touch the front- 
viewing-window. Once completed, the O-ring compression plate was tightened. This 
compressed the O-rings (which were slightly undersized because the tubes were so close 
together) providing a seal between the tube support block and the stainless-steel backing 
plate. Also, the compression plate had grooves for the O-rings to fit into to ensure proper 
alignment and seal. | 

2. System Leakage Check 

After the system was isolated from the atmosphere and system integrity 
restored, a Seargent Welch 10 SCFM vacuum pump was turned on and the pressure in the 
system was reduced to 25 in-Hg vacuum. The system was then secured and left untouched 
for at least 10 hours to see if there was any leakage. If there was significant leakage(>1 
inHg), then the system vacuum was broken by cracking valve R-2 (Figure 2) slowly (to 


ensure minimum moisture was allowed to enter the system). This valve was left open until 
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the system was at atmospheric conditions. The system was then pressurized with air to 
approximately 5 psig through valve R-2. Large leaks were then detected by simply 
listening to the air issuing from the system. If the leaks were small, they could be found by 
applying a small amount of soap/water solution to all joints in the system. One leak that 
was found was where the sleeve had been soldered onto the Turbo-B tube. Since this was 
found to be the source of a leak, all the Turbo-B tube joints were carefully resoldered. 
After all leaks had been corrected, the system was again subjected to a vacuum for a 
minimum of 10 hours. If the increase was less than 1/2 to 1 inHg rise in a 10 hour period, 


then the system was filled with refrigerant. 


E. REFRIGERANT 
1. Fill 
a. From System Storage Tank 

A refrigerant storage tank was used to store excess R-113 during 
experiments. The storage tank prevented discharge of the R-113 into the atmosphere and 
made experimentation less costly. To fill the evaporator with R-113 from the storage tank, 
the system pressure was equalized to the storage container pressure by opening valves R-1 
and R-7 (Figure 2). Once the pressure was equal, then valves R-6 and R-4 were opened. 
Due to the pressure difference caused by the hydrostatic head of the liquid in the storage 
container, R-113 flowed from the storage tank into the evaporator. The amount of 
refrigerant that was transferred was controlled by throttling valve R-4 to obtain the desired 
level. Once the desired amount of refrigerant had been transferred, valves R-4, R-6, R-7, 
R-8 were closed. If required, additional R-113 was later added from a storage drum at 


atmospheric pressure to the system using valve R-2 (see next section). 
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b. From Refrigerant Storage Canister 
To fill the apparatus from the refrigerant drum, the ethylene glycol/water 
mixture was first cooled to a temperature of 5°C. Both condenser pumps were put into 
operation thus allowing coolant to flow through the condenser section. This allowed for a 
differential pressure to exist between the evaporator and ambient. A clean plastic hose was 
connected to valve R-2 (Figure 2) and the other end was end put in the refrigerant drum. 
Once this differential pressure and hose hookup was completed, refrigerant transfer 
between ambient and the evaporator was possible by slowly opening valve R-2. 
2. Removal 
a. To System Storage Tank 
For tube replacement, system maintenance or clean-up purposes, the R- 
113 was transferred to the storage tank. The ethylene glycol/water mixture temperature 
flowing through the storage tank was cooled to less than 5°C; valves R-7 and R-8 (Figure 
2) were opened and the vacuum pump was turned on to put the storage tank under vacuum. 
Once the storage tank was under a vacuum, R-113 flowed from the evaporator to the 
Storage tank because of the pressure difference. This pressure difference was maintained 
by keeping the heaters on in the evaporator and the coolant flow through the storage tank. 
As the level in the evaporator came down, power to the heaters was secured to avoid 
overheating. Any remaining R-113 that could not be transferred into the storage tank was 
drained through valve R-5 and put into a drum to be used later. 
b. To Refrigerant Storage Canister 
This procedure is primarily used to change refrigerants. The system 
vacuum can be broken by slowly cracking open R-2 (Figure 2) to ensure that moisture was 


not allowed to enter the system. Once the system was at atmospheric pressure, valve R-2 
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was opened wide to allow air to enter when the refrigerant was removed. The refrigerant 


could then be drained through valve R-5 and put into a drum to be used later. 


F. OPERATION 
1. System Lightoff, Securing and Emergency Procedures 
See Appendix E 
2. Normal Operation 

The evaporator was filled with R-113 to a level of approximately 10 cm above 
the top tubes in the bundle. Prior to operating the system, the 8-ton refrigeration unit was 
run for approximately an hour to reduce the ethylene-glycol/water mixture in the sump to a 
temperature of 5°C. The pressure in the evaporator was usually 14 to 16 inHg vacuum if 
the system was left secured overnight. As the sump was brought to temperature, the data 
acquisition system and computer were turned on. This allowed the temperature in the 
system to be monitored during warm-up to saturation conditions. With this and pump 
number one running, the auxiliary heaters were then used to bring the pool up to saturation 
temperature. The auxiliary heaters were set to 1 kW nominal after pool temperature had 
been reached. The system pressure was slowly raised by allowing minimal coolant to flow 
at a low rate through the condenser tubes while heating the system with the auxiliary 
heaters. Once the required saturation temperature of 47.5°C was reached, the instrumented 
tubes were switched on and set to the desired value. This was done to prevent the tubes 
from nucleating prematurely. The heat flux of the instrumented tubes was then slowly 
increased by adjusting the rheostat. For increasing heat flux, the data were taken with 
very small] heat flux increments, waiting at least 5 minutes to attain steady state conditions 
at each heat flux. At all regions of the boiling curve (and especially near the onset of 
nucleate boiling), two readings were taken for each heat flux. The bundle was 


continuously monitored through the observation windows. 
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Figure 12 shows the tube bundle arrangements used during experimentation. 
Test One was with only instrumented tube number one on in the bundle. Test two was 
with instrumented tubes one and two on in the bundle. Test three was with instrumented 
tubes one, two and three on in the bundle. Test four was with instrumented tubes one 
through four on in the bundle. Test five was with all five instrumented tubes on in the 
bundle. Test six was all five instrumented tubes plus all five pairs of active heater tubes on 


in the bundle. Test seven was the same as test six but with the simulation heaters on. 


G. DATA REDUCTION PROCEDURES 

The data reduction program "DRP4A" was used during the experiments for 
processing the data collected. The program was written in HP Basic 3.01 and run on an 
HP-9000 series computer. The characteristics and capabilities of this software and the 
entire listing of the program is similar to that provided by Anderson [Ref. 7]. The 
following modifications were made: 

1. Correction for pool height 

2. Installation of new thermocouple at bottom of liquid pool (bundle inlet 
temperature) 

The entire program is included in Appendix F to provide continuity for future thesis 


work. 
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Figure 11. Photograph of Tube Bundle Support Block 
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Figure 12. Tube Bundle Arrangements used During Experimentation 
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V. RESULTS AND DISCUSSION 


A. GENERAL COMMENTS 

The results are presented in five sections with sub-sections as appropriate. The first 
section discusses the natural convection, nucleate pool boiling phenomena and hysteresis 
effects within a small enhanced tube bundle. The second section discusses how evaporator 
pool height above the tube bundle affects the nucleation site activation and general 
evaporator performance. The third section looks at the effects of auxiliary heaters. The 
fourth section discusses the effect of surface history on nucleation site 
activation/deactivation. The final section shows comparisons of data taken during this 
thesis with previously obtained data at the Naval Postgraduate School. 

A list of data runs conducted during this investigation may be found in Appendix A. 
All data files used in this thesis use the following filename sequence. Each file is 


composed of five sets of alpha-numeric characters used to describe the experiment. 


First set (1 char.) Heat flux Increasing heat flux) 
D(Decreasing heat flux) 
Second set (2 char.) Type of tube used TB (Turbo-B) 
Third set (2 char.) Number of Experiment 
for accounting purposes 


To give an example, the file name ITB12 means "increasing heat flux, using a Turbo- 
B tube and the 12th data run. If more detail is desired about a data run, see Appendix A. 
All plot files are similar to data files except they start with the letter "P". All tests were 
done with the auxiliary heaters at either 0, 1, or 3 kW nominal power setting (for details see 


Appendix A). 
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All the graphs show heat flux (W/m2) along the ordinate (y axis) and wall superheat 
(K) along the abscissa (x axis). The heat flux was changed through the use of a variac 
from 600 to 100,000 W/m? for increasing or 100,000 to 600 W/m2 for decreasing runs. 
Approximately 25-30 points were taken on the increasing heat flux and about 20 points 
were taken for the decreasing heat flux. The wall superheat is the difference between the 
corrected average tube wall temperature and the liquid saturation temperature. 

To clarify the operation of the evaporator during the above tests, we used seven basic 
bundle arrangements during the test procedures as discussed and shown in Chapter IV. 

Anderson [Ref. 7] used four types of surface preparation techniques. 

1. Surface preparation A maintains the evaporator tube(s) at a heat flux of 30 kW/m2 
for one hour, slowly reducing heat flux to a minimum value of around 1 kW/m2, and 
beginning data collection immediately at successively increasing heat fluxes. The 
motivation for this surface preparation is a widely varying dynamic load often encountered 
in an air-conditioning machine. 

2. Surface preparation B simulates an operating refrigeration system undergoing 
frequent off-cycles. A heat flux of 30 kW/m2 at saturation conditions is maintained for one 
hour. Power is then secured for 30 minutes and then data are taken at successively 
increasing heat flux steps. 

3. Surface preparation C represents the initial start-up of a refrigeration system. 
During this process, the evaporator was brought up to temperature and pressure with the 
auxiliary heaters alone to ensure the tubes did not have any active nucleation sites. Once at 
Saturation pressure, the tube was turned on and the heat flux varied from 600 to 100,000 
W/m? to develop a heat flux versus wall superheat plot. This process is very interesting as 


it introduces many complexities such as boiling hysteresis. 
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4. Surface preparation D represents a continuous operating air-conditioning system. 
It is closely modelled by a continuously decreasing heat flux starting from maximum power 
of 100 kW/m¢ after it has been boiled at this power for one-half hour. Considering that the 
return flow to the boiler is normally two-phase, nucleation sites are maintained active 
throughout the evaporator tube bundle. The decreasing heat flux data runs are normally the 
most readily reproducible and are used most often in measurements of tube-bundle 
performance. This surface preparation is the most widely used for reporting the 
performance of evaporator-tube surfaces. 

Three of these surface preparation techniques were used during the course of the 
present work. Surface preparation technique C and D were used exactly as stated above 
and surface preparation technique B was the same but with a modification. Surface 
preparation technique B was modified such that power to the tube heaters could be secured 
for any amount of time but the pool temperature had to be maintained at saturation 
temperature. 

1. Repeatability of data 

Figure 13 presents similar data taken on two different days for the top tube 
operating alone in the bundle; the repeatability is seen to be very good. The arrows on the 
figure indicate the familiar hysteresis loop. The two tests show very similar temperature 
overshoot values (indicated by the honzontal line). The point of incipience (point where 
nucleate boiling first occurs) is also very similar for both tests. 

The heat-transfer coefficient, h, is often used as a measure of tube performance: 
it is found by dividing the heat flux by the wall superheat. From Figure 13, values of h in 
the nucleate boiling region range from 1540 W/m2K at low heat flux to 13,750 W/m7?K at 
high heat flux. Corresponding values of h in the convection region range from 480 


W/m2K at low heat flux to 700 W/m2K at medium heat flux. It can be seen from Figure 
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13 that the tube performance typically improves by a factor of five once the tube starts to 
nucleate. This explains the importance of ensuring a flooded evaporator operates in the 
nucleate boiling rather than the natural convection region. 

Also shown on Figure 13 are typical uncertainty bands for the data at both high 
and low heat flux and wall superheat. It can be seen that the greatest uncertainty occurs at 
low wall superheat and heat flux; this was as expected since the values of wall superheat 
are of the order accuracy for the thermocouples. 

Figure 14 shows the limited predictions that are available for both natural 
convection and nucleate pool boiling data. These correlations were formulated for single 
smooth tubes, but represents all that is available for enhanced surfaces and tube bundles. 
This indicates the urgent need for further work in this area. For the natural convection 
region, two correlations are shown: Churchill and Chu [Ref. 24] and Churchill and Usagi 
[Ref. 25] . The experimental data for a single tube within the bundle show an appreciable 
higher heat-transfer performance than prediction. This might be expected as the tube is an 
enhanced surface. However, in the convection region, One should not expect a large 
difference in performance between different tube surfaces as they are designed primarily for 
the nucleate region. Indeed, experiments by Murphy [Ref. 6] have shown that the data for 
enhanced surfaces falls on top of smooth tube data in the convection region. 
Consequently, the discrepancy seen is more likely to be due to circulation effects within the 
bundle. 

In the nucleate boiling region, there is also a large enhancement seen in the 
experimental data when compared with the correlation of Stephan and Abdelsalam [Ref. 
26]. However, in this region, the discrepancy can be explained due to the enhanced nature 
of the surface. Circulation effects are still likely to be present, but small for a single tube 


when compared to the large number of active nucleation sites. 


40 


B. BUNDLE TRENDS FOR A POOL HEIGHT OF 10 CM 

Figure 15 shows the influence of tube two on tube one. One may expect the thermal 
currents from the lower tube to further enhance the upper tube. However, with this lower 
tube heated, the performance of tube one in the convection region did not change. More 
striking is the significantly poorer performance of tube two in this region for which no 
explanation can be found. It is interesting to note that tube two (which should behave as 
an isolated single tube in the convection region, 1.e. feeling no effect from the tube above) 
is now in much better agreement with the natural convection correlations. 

Also seen is that the top tube nucleates at a lower heat flux. However, once both 
tubes are nucleating, the agreement is very close, with the top tube having the slightly 
better performance. Although this slight discrepancy could be due to uncertainty in the 
data, it is more likely to be due to the influence of the bubbles coming from the tube below, 
especially when one compares it with the top tube acting alone (Figure 13). 

Figure 16 shows decreasing heat flux for two tubes. Agreement is very close, with 
the top tube performing slightly better at high heat flux, but the opposite at low heat flux. 
However, as heat flux is reduced, the uncertainty in the data increases. 

With three tubes in operation (Figure 17) , the trends are much the same. One point 
of interest is the effect on lower tubes as upper tubes start to nucleate. As tube one starts 
to nucleate, tubes two and three shift significantly to the left, indicating an increase in 
heat-transfer coefficient. However, tubes two and three remain in the convection region. 
Any further movement to the left by the upper tubes causes a corresponding shift to the left 
for the tubes below. This 1s probably due to increased circulation effects (i.e. as the 
bubbles from the upper tubes move around the evaporator, they help to agitate the thermal 
boundary layer around the tubes below). The shift to nucleation is strictly in order as one 


moves down the bundle. 
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For decreasing heat flux (Figure 18), there is a significant drop-off in performance 
for tube three, especially at a low heat flux. It may be that as one moves down the bundle, 
the greater pressure deactivates the nucleation sites at a faster rate. 

Figures 19 to 22 show four and five tubes activated for increasing and decreasing 
heat flux. Similar trends to those already reported are seen; tube one nucleates first 
followed by two, three, etc. As each successive tube nucleates, the performance of the 
remaining lower tubes is improved slightly. In the nucleate region at high heat flux, tube 
one exhibits the best performance followed by tubes two, four, five and three. This 
somewhat strange order maybe due to uncertainty or the fact that tube three 1s only partially 
nucleating. At the very highest heat flux, all five tubes seem to coincide. At low heat 
fluxes, tube two gives the best heat transfer; it is unclear as to why this should be. 

Because of the unexpected behavior of tube three at high heat flux, the positions of 
tubes three and five were switched to see if this behavior was a "bundle" effect or simply 
due to the particular tube itself. Figure 23 shows the results after switching for decreasing 
heat flux. New tube five (formerly tube three) now has the poorest performance, indicating 
that it is the tube rather than the bundle causing the earlier effect. It may be that air pockets 
surround the thermocouple junctions of tube five. Consequently, the tubes were kept in 
this position for the remainder of the tests, with the "bad" tube in the lowest position. 

Figures 24 and 25 show the entire bundle activated for both increasing and decreasing 
heat flux while Figures 26 and 27 show the entire bundle plus simulation heaters activated 
for increasing and decreasing heat flux. These simulation heaters give the bottom of the 
bundle an inlet quality of about 20 percent, thereby simulating a real evaporator. Similar 
trends to those seen before are evident. The heat flux at which incipience occurs continues 
to decrease for all tubes, indicating the extra quality moving through the bundle helps to 


trigger “early” nucleation. As this simulates part of a real evaporator, this trend is 
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encouraging, as it is more likely to be operating in the nucleate region rather than the 
convection region. For the decreasing heat flux run, the heat transfer from the lower tubes 
has also been improved due to the increased vapor passing over them. 

As already mentioned, it is difficult to make too many direct comparisons when all 
five tubes are independently instrumented. One of the tubes may have certain induced 
errors due to manufacture that another tube does not. Consequently, of perhaps more use 
is to look at what happens to the top tube only as successive tubes are activated below. 
Figure 28 does this for all seven tests (see Chapter IV: normal operation, for description of 
tests) for an increasing heat flux. In the natural convection region, the point of incipience 
occurs at increasingly lower heat fluxes as more tubes are activated. In the nucleate boiling 
region, the heat-transfer coefficient increases as more tubes are activated. Both the natural 
convection and nucleate boiling regions are well defined indicating good reliability of the 
data. For decreasing heat flux (Figure 29), the heat-transfer coefficient 1s again clearly 
increased as more tubes are activated below. This is more evident at low heat flux than at 


high heat flux. 


C. BUNDLE PERFORMANCE WITH VARYING POOL HEIGHT 

One of the objectives of this thesis was to see the effect of pool height on heat- 
transfer performance. Pool height here is defined as the height from the free liquid surface 
down to the top of the bundle (tube one). All experiments to date had been carried out at a 
pool height of 10 cm. Two additional heights of 0 cm and 20 cm were also studied with 
the same vapor pressure within the evaporator (i.e. to see the effect of pressure head on the 
tubes in the bundle). It is understood that a true flooded evaporator has zero pool height 
(indeed the top tubes operate in the bubbly froth). 

Figures 30 to 42 show the same seven tests for a pool height of 20 cm as previously 


shown for 10 cm. The trends of all these graphs are exactly the same as for 10 cm. The 
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effect of the greater pool height is generally to delay the point of incipience. Figures 43 and 
44 show tube one only for all seven tests (i.e. to see the effect of successive activated tubes 
below) for increasing and decreasing heat flux respectively. This to shows similar trends 
to before, with the effect of successive tubes triggering earlier nucleation on the top tube as 
well as enhancing the nucleate boiling heat-transfer coefficient. 

Figures 45 to 57 show the same seven tests for a 0 cm pool height and again the 
trends are very similar. However, this time, the point of incipience occurs at the lowest 
value of heat flux. Figures 58 and 59 show tube one only for all seven tests for increasing 
and decreasing heat flux respectively. The conclusions drawn from these figures are the 
same as before. 

To summarize the effect of pool height, Figures 60 and 61 show tube one (i.e. Test 
One ) at the three pool heights for increasing and decreasing heat flux. It can be seen that 
the zero pool height data nucleates first followed by 10 cm and then 20 cm. Similar trends 
were seen regardless of what test was chosen 1.e. for test seven with whole bundle plus 
simulation heaters activated, all three points of incipience (for the 3 pool heights) move to a 
lower heat flux but the zero pool height data still nucleate first. Obviously, the greater 
pressure heads felt by the tubes at larger pool heights inhibit the transition to nucleate 
boiling 1.e. larger wall superheats are required to overcome the increased pressure. Further 
experiments should be conducted in the future whereby the vapor pressure above the 
bundle is varied with the pool height such that the pressure that each tube "feels" at a given 


location is the same. 


D. EFFECT OF AUXILIARY HEATERS 
For all tests reported so far, the auxiliary heaters were set to 1kW nominal. The 
auxiliary heaters were needed to maintain pool temperature throughout the test. They were 


only really needed at low heat fluxes or if only a few tubes in the bundle were activated. 
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For these tests, if they were not used, the evaporator could not produce enough vapor to 
keep the pressure constant even with the condenser shut down (i.e. the walls of the test 
section acted as the condenser). As a consequence, the pool temperature would drop. 
However, the necessary operation of the auxiliary heaters has an effect on the circulation 
patterns within the bundle (this has been demonstrated in the single tube facility). 
Therefore, by having these auxiliary heaters set to the same value for each test, any effect 
they may have should be a constant factor. To see the effect that the auxiliary heaters have, 
Figure 62 shows data for three tubes on during a decreasing heat flux with the auxiliary 
heaters off compared to Figure 18 with them set at 1 kW for a 10 cm pool height. 
Essentially the effect 1s minimal in the nucleate region at high heat fluxes. At lower heat 
fluxes, no data could be taken with the auxiliary heaters off. In the convection region 
(again at low heat flux so there is no data with the auxiliary heaters off), the effects of 
increasing the auxiliary heater power is to increase the heat-transfer coefficient (i.e. 
reducing the wall superheat). The same effect was seen for four tubes with the auxiliary 
heaters off and set to 1 kW for a 10 cm pool height; shown in Figures 63 and 20 
respectively. Figures 64 and 30 show the same trend for tube one with the auxiliary 
heaters set at 3 kW compared to 1 kW for a pool height of 20 cm. These trends are 
probably seen due to the increased circulation within the bundle at the higher auxiliary 
heater setting. Obviously, it would be preferable to do away with the auxiliary heaters 


altogether. With R-114, due to the lower normal boiling point, this should be possible. 


E. NUCLEATION SITE ACTIVATION/DEACTIVATION 
Nucleation site activation/deactivation is not well understood. Once a tube is 
nucleating, there are a couple of ways that the nucleation sites can die out: 
1. Reducing the heat to the tube surface thereby allowing the fluid to collapse the 


vapor bubble and re-wet the surface due to pressure effects. 
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2. Increasing the subcooling of the liquid around the tube (i.e. causing the vapor 
bubble to condense into the fluid). 

The following tests were conducted to see how quickly the nucleation sites, once 
activated, died away if the power to the tube was cut off (i.e. to see how quickly the liquid 
re-wets the surface as discussed in mechanism one above). The pool temperature was 
maintained at saturation conditions throughout the test (i.e. there was no subcooling). All 
the tests were conducted on the top tube only. The test procedure used was the modified 
surface preparation B. The power to the heater was then turned off for a known finite time 
and then reactivated. An increasing run was then carried out. This was repeated for a 
number of different times and pool heights. The auxiliary power setting was 1 kW for all 
tests. 

Figures 65 and 66 show the data for a pool height of 20 cm. The turn-off time for the 
heater varied between five seconds and one hour. It can be seen that after only 15 s, the 
nucleation sites have all died out and the tube behaves as if it had been secured overnight 
(these are the data sets listed as 14 and 15 hours). Of course, by securing overnight, 
subcooling also plays a part in deactivating the sites. 

Figure 67 shows the data for a pool height of 10 cm. The turn-off time for the heater 
varied between five seconds and one minute. It can be seen that after one minute, the 
nucleation sites have all died out and the tube behaves as if it had been secured overnight. 
Hence, it would seem that due to the reduced pressure head trying to push the liquid into 
the sites, the time needed to deactivate all the sites is significantly longer. 

Figures 68 and 69 show the data for a pool height of 0 cm. The turn-off time for the 
heater varied between 40 seconds and 8 hours (keeping pool temperature constant 
throughout). It can be seen that even for a shut-off time of 8 hours, none of the nucleation 


Sites are deactivated. The pressure head is negligible here and it is mechanism two that 
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“deactivates the sites. As the pool temperature was not subcooled during this particular test, 
none of the sites are therefore deactivated. Figure 70 shows a test where the apparatus was 
secured for 2.75 hours (i.e. allowed to subcool). The deactivation seen is due to the 
subcooling that has occurred. 

Similar trends for all three pool heights have been obtained for an auxiliary heater 
power setting of 3 kW. Slight differences in the times needed for total deactivation were 
noted to be 20 seconds for 10 cm vice one minute. These data are shown in Figures 71 to 


fiey 


F. COMPARISONS WITH PREVIOUS NPS DATA 

Figure 76 shows a comparison between the present data for a Turbo-B bundle in R- 
113 and the data of Anderson [Ref. 7] for a smooth tube bundle in R-113. For clarity, 
only the top tube has been shown for a decreasing heat flux. It can be seen that the heat 
transfer enhancement given by the enhanced tube is three to five times for high and low 
heat fluxes respectively. 

This is also true for the overall bundle heat transfer coefficient, giving a factor of 4 
enhancement at a typical operating heat flux of 20,000 W/m2. Also shown on the figure is 
the correlation of Stephan and Abdelsalam [Ref. 26] showing a fairly good agreement with 
the smooth tube data. 

Figure 77 shows a comparison between the present data and that of Eraydin [Ref. 9] 
who used a Turbo-B bundle in R-114. It can be seen that there is very little difference 
between the two sets of data. The R-114 data giving the slightly better performance. This 
could be due to the slightly better thermal conductivity of R-114 at its saturation 


temperature compared to R-113 at its saturation temperature. 
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Figure 13. Performance of Tube 1 at a 10 cm Pool Height for 
Similar Data Taken on Different Days to Show Repeatability 
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Figure 14. Performance of Tube 1 at a 10 cm Pool Height Along With 
Correlations for Natural Convection and Nucleate Pool Boiling Regions 
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Figure 15. Performance Variation of Tubes 1&2 at a 10 cm Pool 


Height During a Increasing Heat Flux (ITB05) 
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Figure 16. Performance Variation of Tubes 1&2 at a 10 cm Pool 


Height During a Decreasing Heat Flux (DTB05) 
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Figure 17. Performance Variation of Tubes 1 to 3 
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Figure 18. Performance Variation of Tubes 1 to 3 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB06) 
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Figure 19. Performance Variation of Tubes 1 to 4 at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB07) 
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Figure 20. Performance Variation of Tubes 1 to 4 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB07) 
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Figure 21. Performance Variation of Tubes 1 to 5 at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB27) 
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Figure 22. Performance Variation of Tubes 1 to 5 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB27) 
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Figure 23. Performance Variation of Tubes 1 to 5 at a 10 cm Pool 
Height During a Decreasing Heat Flux with Tubes 3&5 Switched (DTB30) 
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Figure 24. Performance Variation of the Bundle at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB48) 
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Figure 25. Performance Variation of the Bundle at a 10 cm Pool 


Height During a Decreasing Heat Flux (DTB48) 
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Figure 26. Performance Variation of the Bundle plus Simulation 
Heaters at a 10 cm Pool Height During a Increasing Heat Flux (ITB49) 
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Figure 27. Performance Variation of the Bundle plus Simulation 
Heaters at a 10 cm Pool Height During a Decreasing Heat Flux (DTB49) 
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Figure 28. 10 cm. Performance Variation of Tube 1 in a 
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Turbo-B Bundle During a Increasing Heat Flux 
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Figure 29. 10 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influence by a Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Decreasing Heat Flux 
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Figure 30. Performance of Tube I at a 20 cm Pool Height 
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Figure 31. Performance Variation of Tubes 1&2 at a 20 cm Pool 


Height During a Increasing Heat Flux (ITB34) 
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Figure 32. Performance Variation of Tubes 1&2 at a 20 cm Pool 


Height During a Decreasing Heat Flux (DTB34) 
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Figure 33. Performance Variation of Tubes 1 to 3 at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB35) 
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Figure 34. Performance Variation of Tubes 1 to 3 at a 20 cm Pool 


Height During a Decreasing Heat Flux (DTB35) 
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Figure 35. Performance Variation of Tubes 1 to 4 at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB36) 
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Figure 36. Performance Variation of Tubes 1 to 4 at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB36) 
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Figure 37. Performance Variation of Tubes 1 to 5 at a 20 cm Pool 


Height During a Increasing Heat Flux (ITB37) 
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Figure 38. Performance Variation of Tubes 1 to 5 at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB37) 
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Figure 39. Performance Variation of the Bundle at a 20 cm Pool 


Height During a Increasing Heat Flux (ITB38) 
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Figure 40. Performance Variation of the Bundle at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB38) 
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Figure 41. Performance Variation of the Bundle plus Simulation 
Heaters at a 20 cm Pool Height During a Increasing Heat Flux (ITB40) 
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Figure 42. Performance Variation of the Bundle plus Simulation 
Heaters at a 20 cm Pool Height During a Decreasing Heat Flux (DTB40) 
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Figure 43. 20 cm. Performance Variation of Tube 1 in a P 
ool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 


Turbo-B Bundle During a Increasing Heat Flux 
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Figure 44. 20 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Decreasing Heat Flux. 
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Figure 45. Performance of Tube 1 at a 0 cm Pool Height 
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Figure 46. Performance Variation of Tubes 1&2 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB52) 
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Figure 47. Performance Variation of Tubes 1&2 at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB52) 
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Figure 48. Performance Variation of Tubes 1 to 3 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB53) 
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Figure 49. Performance Variation of Tubes 1 to 3 at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB53) 
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Figure 50. Performance Variation of Tubes 1 to 4 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB55) 
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Figure 51. Performance Variation of Tubes 1 to 4 at a 0 cm Pool 


Height During a Decreasing Heat Flux (DTBSS) 
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Figure 52. Performance Variation of Tubes 1 to 5 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB56) 
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Figure 53. Performance Variation of Tubes 1 to 5 at a 0 cm Pool 


Height During a Decreasing Heat Flux (DTB56) 
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Figure 54. Performance Variation of the Bundle at a 0 cm Pool 


Height During a Increasing Heat Flux (ITB57) 
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Figure 55. Performance Variation of the Bundle at a 0 cm Pool 


Height During a Decreasing Heat Flux (DTB57) 


90 


100000 

















10000 “HHH 
Seer eer MERU Goa onli Gee Ee RECT Oe Pee OE Ae Ses 
cee | ee A ee Dt eee ee a 
= a 2 dd th ob ep 
: eee ee Lee eae 
S Ct mea 
BS | Ta 
ka 
© : 
. | 
1000 Bill 
eer | 
eer 
foo 
as 






AHHH all 


(Twall-Tsat)/(K) 


Figure 56. Performance Variation of the Bundle plus Simulation 
Heaters at a 0 cm Pool Height During a Increasing Heat Flux (ITB58) 
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Figure 57. Performance Variation of the Bundle plus Simulation 
Heaters at a 0 cm Pool Height During a Decreasing Heat Flux (DTB58) 
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Figure 58. 0 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Increasing Heat Flux 
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Figure 59. 0 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Decreasing Heat Flux 
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Figure 60. Comparison of Tube 1 in a Pool of R-113 for Three 
Different Pool Heights During a Increasing Heat Flux 
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Figure 61. Comparison of Tube 1 in a Pool of R-113 for Three 
Different Pool Heights During a Decreasing Heat Flux 
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Figure 62. Performance Variation of Tubes 1 to 3 at a 10 cm 
Pool Height During a Decreasing Heat Flux with Auxiliary 
Heaters off (DTBO6T) 
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Figure 63. Performance Variation of Tubes 1 to 4 at a 10 cm 


Pool Height During a Decreasing Heat Flux with Auxiliary 
Heaters off (DTBO07T) 
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Figure 64. Performance of Tube 1 at a 20 cm Pool Height with the 
Auxiliary Heater Setting at 3kW 
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Figure 65. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 20 cm Pool Height 
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Figure 66. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 20 cm Pool Height 
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Figure 67. Nucleation Site Deactivation with a 1 kW Auxiliary 


Heater Setting and a 10 cm Pool Height 
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Figure 68. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 0 cm Pool Height 
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Figure 69. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 0 cm Pool Height 
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Figure 70. Nucleation Site Deactivation with a 1 kW Auxiliary 


Heater Setting and a 0 cm Pool Height 
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Figure 71. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater Setting and a 10 cm Pool Height 
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Figure 72. Nucleation Site Deactivation with a 3 kW Auxillary 
Heater Setting and a 10 cm Pool Height 
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Figure 73. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater setting and a 10 cm Pool Height 
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Figure 74. Nucleation Site Deactivation with a 3 kW Auxiliary 


Heater Setting and a 10 cm Pool Height 


109 


100000 - — ee 


10000 


Heat Flux/(W/m*2) 


1000 











O% Oi iit 
Turbo-B Bundle| | || Fh 





100 


Twall-Tsat/( K) 


Figure 75. Nucleation Site Deactivation with a 3 KW Auxiliary 


Heater Setting and a 10 cm Pool Height 
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Figure 76. Comparison of Tube 1 During a Decreasing Heat Flux for 
Smooth vs Turbo-B Bundles in R-113 
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Figure 77. Comparison of Tube 1 During a Decreasing Heat Flux for 
Turbo-B Tubes in a Bundle with R-113 vs a Bundle for R-114 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

1. When secured overnight, the bundle exhibited significant temperature overshoot 
for increasing heat flux (approx. 10°C). The hysteresis loop formed by increasing 
followed by decreasing heat flux was very repeatable. 

2. In the natural convection region, activating lower tubes within the bundle did not 
affect the upper tubes. However, the heat-transfer coefficient decreased as one moved 
down the bundle. 

3. Activation of the lower tubes caused the upper tubes to nucleate earlier (i.e. at a 
lower heat flux). As each tube successively nucleated, the heat-transfer performance of the 
remaining lower tubes improved. 

4. In the nucleate boiling region, the performance of the top tube was the best. This 
performance dropped off steadily as one moved down the bundle. 

5. Variation in liquid height above the bundle (from 0 to 20 cm) significantly delayed 
the point of incipience. 

6. For zero pool subcooling, nucleation site deactivation was greatly dependent on 
both pool height (i.e. pressure head) and tube turn off tme. At zero pool height, 
deactivation was negligible for all ttme whereas for 10 and 20 cm, there was complete 
deactivation in 60 and 15 seconds respectively (i.e. decreasing with increasing pool 
height). 

7. Ata typical heat flux of 20 kW/m2, the Turbo-B tube bundle heat-transfer 


coefficient was approximately four times that of a smooth tube bundle. 
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8. The uncertainty of the data was greatest at low heat flux and low wall superheat 


due to the measurement being of the order of accuracy of the thermocouples. 


B. RECOMMENDATIONS FOR FUTURE WORK 

1. Due to the way a real flooded evaporator operates, it would be preferable to heat 
the tubes within the bundle with hot water. 

2. The pool height tests should be repeated, but by varying the vapor pressure above 
the liquid, the local pressure at the tube surface should be maintained at a constant value. 

3. A pressure test should be conducted and a pressure relief valve fitted to enable the 
testing of alternative refrigerants such as R-124 which has a higher vapor pressure. 

4. Attention needs to be given to the question of refrigerant disposal. There are 
reclamation projects undertaken by most manufacturers; however, a method still needs to 
be found to remove the refrigerant from the apparatus into a container for reclamation. 

5. Neutrally buoyant particles should be placed in the pool to facilitate study of the 
circulation patterns within the bundle in more detail. 

6. A high speed camera should be used to study the nucleation process and 
circulation patterns in more detail (the mechanical engineering department is hoping to 
purchase such a camera in the near future). 


7. A sight glass should be fitted to the refrigerant storage tank. 
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10. 


11. 


t2: 


13. 
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APPENDIX C: SAMPLE CALCULATIONS 


Data set number 1 of run DTBO3 (Decreasing heat flux, Turbo-B tube, Surface 
preparation D) was used for the sample calculations in order to validate the program used 


for data acquisition DRP4A. The working fluid was R-113. 


1. Test tube dimensions 
Dc = 11.60 mm 


D, = 14.15 mm 
D; = 12.70 mm 
L =203.2 mm 
Lu = 25.4 mm 


2. Measured Parameters 
T1 =54.41 °C 
T2 = 55.25 °C 
T3 = 55.80 °C 
T4 =55.15 °C 
TS = 55.04°C 
T6 = 55.38°C 
Tld1 = 47.46°C 
Tld2 = 47.48°C 
Aas = 3.838V 
Vas = 3.958V 
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3. Calculations 


The heater power is first calculated from 
q=Vas(V) x Aas(V) x 60(V/V) x 1(A/V) 


Note: The multiplication factors of volt and amp sensors are 60 and 1, respectively. 


Therefore: 
q= 3.958V(3.838V)(60V/V)(1A/V) 
= 911.45W 


The tube inside wall temperature is obtained from the average of all six thermocouple 


readings. 


= 1/6 (54.41 + 55.25 + 55.80 + 55.15 + 55.04 + 55.38) 
= 55.17°C 


The tube outside temperature is calculated by knowing the inside wall temperature and 
using Fouriers Conduction law . Uniform radial conduction is assumed. 


q inS*)) 


271 Key)(L) 


Two = Twi - 
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where the second term on the right hand side is the Fourier conduction term. If we 


define this term as 
D 
In(—2 
ae ql 5 
2M key)(L) 


and 


6, = Two - Tsat, 


where kcu is the thermal conductivity and is calculated as follow 
kcu = 434.0 - [0.1(Twi-K)] 
= 434.0 - [0.1(328.32)] 
= 401.17 W/mK 


Now, 
14.15 
911.45[InG on 7 0) 


27(401.17)(0.2032) 
= (55.17 - 0.3567)°C 
= 54.82°C 


Two = Twi - 


The liquid saturation temperature at the top of the tube bundle is 


Tsat = ene 


— 47,.46+47,48 
2 


= 47.47°C 
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In order to calculate the local saturation temperature for each tube, correction factors 
are needed to account for hydrostatic pressure differences between the tube locations and 


the liquid free surface. This difference is calculated by: 


AP = p (g)(ht) 
where the density is calculated by 


Pp = 1620.7479 - Tsat(2.2186346 + Tsat(2.3578291E-3)) 


Therefore, for Tsat = 47.47°C 


p = 1510.1 kg/m3 


So, for the top tube in the bundle which is 0.0124 m below the thermocouples 


measuring the pool temperature, 


AP = 1510.1 (9.81)(0.0124) 
AP = 183.69 Pa 


For a 183.69 Pa pressure difference,the corrected saturation temperature is obtained 
by adding 0.054°C (obtained from standard tables for R-113 or by equations shown in 
Appendix F.) to Tsat. The corrected Tsat is: 


Tsatc = (Tsat + 0.054)°C 
= (47.47 + 0.054)°C 


Tsate = 47.52°C 
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Therefore, the wall superheat can be obtained by the following 


6, = Two - Tsat, 
= (54.82 - 47.52)°C 
= 7.30°C 


Now that the wall superheat is known, we need to calculate the heat flux and the heat 
transfer coefficient. To do this, we know that the tube is 12 inches long and is heated in a 
eight inch center portion of the tube. The unheated lengths of the tube are a one inch and a 
three inch section on opposite ends of the tube. These unheated lengths have a fin effect 
during the heat transfer process to the evaporating refrigerant. In order to account for this, 
the following procedure was adopted for both the one and three inch sections. Calculations 
are shown below for the one inch section . 

Heat transfer from the unheated end is calculated as heat from the base of a fin: 


qr = (hy Xp X key X Ac)” X Op X tanh(n x Lg) 


where 
p=1zxD, 
= 1x 14.15E-3 m 
= 44.45E-3 m 
now 
Ac =F (DB - Di) 


Ac = (0.014152 - 0.012707) 


A. = 3.0578E-5 m? 
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The corrected length of unenhanced surface at the end was calculated as follows 


Le=ly +i 
L. = 0.0254 + owe 
L, = 0.0258 m 


hp is the natural convection heat transfer coefficient of the fin like ends and was 
calculated by using the Churchill-Chu [Ref. 24] correlation for natural convection from a 


smooth horizontal cylinder, as modified by Pulido [Ref. 27]. 


Z 


] 
fe ah Secale hf 
L 
hp = 1 0.6 +0.387 x —— ; 
O 


1+ (0.559) 6]27 





L 


where 


| mescr | 0.5 
Koy X Ac 


Therefore an iterative technique was necessary to calculate hp. The iterative technique 
used was to assume hp was 190 ,,2¢ and to continue the iteration until successive values 


are within 0.001 of each other. The fluid physical properties are calculated at the vapor 


mean film temperature, given by the following equation. 


Nee, 


Tsat. + T 
Tram = psa + Two 


Tra = 42:52 + 54.82 


Thm = 51.17°C =324.32°K 
For R-113, the physical properties are given in the program by 


Dynamic viscosity, Tri, in °K 
= 1.34E-5 x Me Eow ess 
Lb = 4.88E-4 kg/m:s 
Specific heat, T,),, in °C 
Cp = 929.0 + (1.03 x Trim) 
Cp = 981.71 J/kg-K 
Density, Tfiim in °C 
P = 1620.7479 - Trim 2.2186346 + T5(0.0023578291)] 
p = 1501.1 kg/m 
Thermal conductivity of R-113, T,),, in °C 
k = 0.071 - (0.000261 x Tgim) 
k = 5.765E-2 W/mK 


Prandtl number 
= ty 
Pr = CG x 7 
= 4.88E-4 
lil la 595; 
Pr = 8.31 
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Thermal expansion 
B = -(1/p)(Ap/AT) 
51.07 = 1501.29 kg/m? 051.27 = 1500.8 kg/m3 
B = -(1/1501.1)(0.49/0.2) 
B = 1.63E-3 (1/K) 


Kinematic viscosity 


Vv = 3.25E-7 m2/s 
Thermal diffusivity 


Qa= 
px Cp 





oy = __5:7645E:-2 
1501.1 x 981.71 
a = 3.912E-8 m2/s 


Knowing the above properties, the heat transfer coefficient, hp, can be obtained by 
iteration. 


hp = 185.42 W/m2-K 


knowing this we calculate n 


r-| be | 0.5 
key XAc 
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185.42 x 44.45E-3 | °° 
401.17 x 3.0578E-5 


N= 25,92 


then we can obtain q¢ 
qr = (hy Xp X Key X Ac)” X Op X tanh(n x Le) 
qr = (185.42 x 44.45E-3 x 401.17 x 3.0578E-5)°> x 7.3 x tanh (25.92 x 0.0258) 
q¢ = 1.36 W 


The corresponding results for the three inch section are 
hp = 158.27 W/m2:K 
dr = 2.04 W 


Therefore the heat transfer through the heated length of the tube is 
ds=q- qf(1 inch section) - qf(3 inch section) 
Gs = (911.45 - 1.36 - 2.04) W 
ds = 908 W 


and the heat flux and the heat transfer coefficient are as follows 


q’ = qs/As 
= q;/(m x Do XL) 
= 908/(x x 0.01415 x 0.2032) 
= 908/9.033E-3 
= 1,005E+5 W/m2 
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and finally the heat transfer coefficient. 


a qs 
A, X (Two-Tsat) 
= 908 / (9.033E-3 x 7.3) 
h = 1.376E+4 W/m2-K 
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APPENDIX D: UNCERTAINTY ANALYSIS 


The same data run (DTBO3) was chosen for the uncertainty analysis. Therefore, the 
measured and calculated parameters found in the sample calculation were used in this 
section. The uncertainty analysis performed was for a high heat flux but the same 
procedure could be performed at any heat flux to determine the uncertainty bands. All 
uncertainties are presented as a percentage of the calculated parameter. The uncertainty 
associated with the experimental parameters is calculated from the equation suggested by 
Kline and McClintock [Ref 28]. For example, if 


R=R (x1, X2,-.-, Xn) 


then 
2 2 2) 0.5 
5R = ae) ee) ‘pene (eas) 
OX) 0x2 OXp 
where: 


OR = uncertainty of the desired dependant variables 
Xn = measured variables 


Ox, = uncertainty in measured variables 


The boiling heat transfer coefficient is given by 


ee 
A, x (Two-Tsat) 
where 
Do 
= _  4q([InG=-)] 
Soe ee Diem 
Two = Twi Pnikcu\(L) 
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In the above equation, the second term on the right hand side is usually called the 


Fourier heat transfer conduction term. If we define this as ®, then 


Do. 
_ q [nts J 
2mkcu)(L) 
and 


@, = Two - Tsat, 


With this notation, the uncertainty in the heat transfer coefficient is obtained using the 


following equation. 




















where 

q=VxI 

q= V(V) x I(V) x 60(V/V) x 1(A/V) 
and the uncertainty is 

8q _[[av]? ,[ar}?]2 

q - [| + |B ? 


The accuracy in the voltage and current sensors are as follows 


dVas = 40.05 V dAas = +0.025 A 


14] 


From the sample calculation section 


Vas = 3.958V Aas = 3.838V 


Therefore, 


5q/q = ((6Vas/V as)” + (SAas/Aas)’) sad 
Sq/q = ((0.05/3.958) + (0.025/3.838)2) °° 
5q/q = 1.42 percent 
Calculation of the surface area and the uncertainty of it are done as follows 
As =%xD, XL 
2 2)+ 
OAs = [a2] an aL | Z 
AS 


Knowing the dimensions of the tube from the manufacturer and estimated 


inaccuracies from work shop tools and human error, the uncertainty was calculated. 


Dimensions 


Do = 14.15 mm L = 203.2 mm 


Inaccuracies in measurements 


5D, =0.1 mm dL = 0.2 mm 
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Uncertainty analysis performed 
BAVAs = ((SDo/Do)" + (SL/L)) 
SAS/Ag = ((0.1/14.15)2 + (0.2/203.2))> 
dA,/A, = 0.7135 percent 


The uncertainty calculation for the Fourier conduction term is given below 
a2 ("fea 
@ q kcu L 
kcu was calculated using 
kcy = 434.0 - [0.1(Twi-K)] 
= 434.0 - [0.1(328.32)]} 
= 401.17 W/mK 
and its uncertainty 


dkcu = [0.1 x 8Twi-K)’] > 


5Twi and 5Tsat are obtained using uncertainties in the thermocouple readings. 
Average wall inside temperature Twi was obtained taking the average of six thermocouple 


readings inside the tube. The uncertainty associated with this variable is 





where 5T, was obtained as + 0.5°C from "Omega 1987 Complete Temperature 
Measurements Handbook pp. T-37. Using engineering judgement this seems high will use 


0.1°C. 
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5Twi = 6 1) ’ ° 
5Twi = 0.04 °C 


Saturation temperature was obtained by taking the average of two thermocouple 


readings and the uncertainty in this temperature was calculated from the following equation. 


ST, ) 
2 
dTsat = 2 x [O.1 ’ 


dTsat = 0.07°C 








dTsat =|2 x 


L 
2 


Knowing the uncertainty in the temperatures, we can now calculate the uncertainties 
in the following. 
Skcu =|(0.1 x 8Twi-K)’| 2 
Skcu =[(0.1 x 273.19)"]2 
dkcu = 27.32 W/mK 


Now we can calculate the uncertainty in the Fourier conduction term 


*aca)* [a]? 





Nd |= 


ie x [0.0142] 2 + | 27.32-)* 10.25! : 
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O® — 6.96 percent 
® 


from the sample calculations we know that ® = 0.3567 °C or 


5 = 0.3567 x 6.96 percent = 0.025 °C 


therefore: 


1 
an» a)” Bs ? [Swi]? , [so], stat"? 
ae be: @ | 16 Op 




















= = [0.0142) 2 + [0.07135] 7 + [o.04 2 [o.02s; a 0.07} oF 


= = 7.4 percent 


Table 3 shows the results of the uncertainty analysis performed. The high and low 
heat flux correspond to the approximate values of 1.005E+5 W/m2 and 597 W/m2 


respectively. 
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TABLE 3. UNCERTAINTY ANALYSIS RESULTS 








Variable igh Heat Fl Low Heat Flux 
AT 7.3 0.38 
Twi 55.17 48.04 
a 47.47 47.6 
éV. 

V 1.26% 16.4% 
ot 

I 0.65% 8.3% 
oq 

q 1.42% 18.% 
5D, 

Do 0.707% 0.707% 
dL 

L 0.098% 0.098% 
5A, 

As 0.714% 0.714% 
OKeu 

kou 6.8% 6.8% 
5D 

@ 6.96% 21.9% 
dh 

h 7.4% 22.4% 
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APPENDIX E: OPERATING PROCEDURES 


YSTEM_ _LIGHT-OF 

1. Power to the 8 ton refrigeration unit is provided by the breakers located in 
the main distribution panel located in the laboratory. These breakers were never secured. 
However, if power to this panel was lost then these breakers must be reset. 

2. Turn the switch on the refrigeration unit control panel, located in front of the 
refrigeration unit to the "auto" position after passing through "on" position. This switch is 
always left on, unless unit was taken down for long repairs. 

2. Push the start button in the control box for the recirculation pump. This 
control box is located on the bulkhead above the recirculation pump in the outside area 
adjacent to the refrigeration unit. 

3. Set the desired temperature on the roughly graduated Fahrenheit scale on the 
control panel thermostat. It requires approximately one hour to chill the sump 5°C. The 
thermometer of the ethylene glycol sump must be monitored to ensure the desired sump 
temperature is attained and maintained. Slight adjustments in the refrigeration unit 
thermostat can be expected due to the coarseness of its scale. 

4. Energize the desired pumps by switching on the breakers in the main 
distribution power panel. Once the power is energized in the main power panel turn on the 
pump motors by pushing down on the arm of the appropriate breaker box for the pumps 
located on the bulkhead next to the ethylene glycol sump. The pumps are marked 
“auxiliary condenser"(Pump #2) and “instrumented tube condenser"(Pump #1), 
respectively in the breaker box. Flow in the auxiliary condensate system can be controlled 
with the individual globe valves located at the coil penetrations on the apparatus. The 


auxiliary condenser will produce the fastest adjustments to system pressure if pressure is 
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rising too quickly. Flow through the instrumented tubes can be controlled by the ball 
valves located at the bottom of the respective flow meter. 

5. Energize the heater variacs desired by switching on the breakers (Auxiliary, 
Simulation, and/or Bundle) in the main distribution power panel (located near entrance of 
lab) and individual breakers for each of these in the power distribution box (near 
apparatus). 

6. After ensuring that the breakers for the heater tubes desired are in the "on" 
position, adjust variac/s to the desired amount of heat flux. Monitor apparatus pressure 
through the pressure gage at the top of the apparatus, ensuring system pressure does not 


exceed 30 psig. 


SYSTEM SHUTDOWN 

1. Turn all variacs to the zero position and switch off all breakers in the power 
panels. 

2. If apparatus will not be operated for an extended period, turn the switch on 
the refrigeration control panel to the "off" position after passing through "on". 

3. Allow the recirculation pump to operate for at least five minutes after 
switching off the refrigeration unit to dissipate any back pressure in the system. 

4. Turn the breakers for the pumps to the off position at the switch boxes near 


the ethylene glycol sump, and then secure power at the main distribution power panel. 


EMERGENCY SHUTDOWN 


1. Secure all power at the main distribution power panel 
2. Pull fire alarm 
3. Evacuate building 
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APPENDIX F: PROGRAM DRP4A 


1000! FILE NAME: ORP4A 

1004! DATE: November 22, 1988 

1008! REVISED: May 1991 (S. MEMORY) 

1@12! 

1@16 G8EEP 

102@ PRINTER IS 1 

1024 Idp=2 

1Q@28! 

1032 PRINT USING “4X,""Select option default is @.""" 
1036 PRINT USING "6X,""@ Taking data or re-processing previous data""" 
104@ PRINT USING "6X,""1 Plotting data on Log-Log 
1044 PRINT USING “6X,""2 Plotting data on Linear""" 
1048 PRINT USING “6X,""3 Purge""" 

1@S2 PRINT USING “6X,""4 Fixup""" 

1@S6 PRINT USING “6X,""S Move""" 

1Q@6Q@ PRINT USING "6X,""6 Comb""" 

1064 PRINT USING "6X,""7 Read Plot""" 

1068! 

1072! IDP IS A PROGRAM VARIABLE TO SELECT A SUBROUTINE 
1076 INPUT Idp 

108@ IF Idp=@ THEN CALL Main 

1084 IF Idp=l THEN CALL Plot 

1088 IF Idp=2 THEN CALL Plin 

1992 1F Ide=3 THEN CALL Purg 

1096 IF Idp=4 THEN CALL Fixup 

110@ IF Idp=S THEN CALL Move 

1104 IF Idp=6 THEN CALL Comb 

1108 IF Idp=7 THEN CALL Readplot 

fitZ2 END 

1116! 

112@ SUB Main 

1124! ICAL=THERMOCOUPLE CALIBRATION 

tiza5 COMS/ Ces CC7) 

1132 DIM Emf(35),7T635),01la(6),02a(6),D1a(6),Doa(6),La(6),Lua(6),Kcua(6),Et(19), 


Ldtc(20),Volt(2),Amp(11),Twa(S),Tw(S), Theta(S), Thetab(S),0(5),Q1(5),Qdp(S) 


VIG 
114Q! 
1144! 
1148 
ise 
1156 
1169! 
1164! 
1168 
fe gi 
176 
118Q! 
1184 
1188 
PVs 
i136 
122 
1204 
1208 
| 


bo 16 


DIM Htube(S) 


THERMOCOUPLE ARRAY (C( )) INITIALIZATION 

Pale O27 1OG8509! 25727. 94569,—-767345.8295, 73025595. 81 
ORM -9247485599,6. 97588641) -2 bel 92E+l Ss, 35. 9407S8Er1 4 
READ C(#) 


PRINT HEADER AND INITIALIZE TIME CLOCK 

PRINTER IS 701 

BEEP 

TINFUT “ENTER MONTH, DATE AND TIME “MM: OD HH-MM 355)" Cares 
OUTPUT DIRECTED TQ DATA AQUISITION SYSTEM ‘HP 3497A) 
OUTPUT 709: “TD";:Dates 

OUTPUT 709: "TD" 

ENTER 709: Dates 


PRINT 
ER INt Month, date and time :“iDates 
PRINT 
PRINT USING “10X,°"NOTE- Program name DRE OQ 7-7 
BEEP 
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1220! DON IS THE VARIABLE FOR DISC NUMBER FOR RECORD KEEPING ONLY 
1224 INPUT “ENTER DISK NUM8ER”,Dn 

1228 PRINT USING “16%, “Uisk numGer Game” (27°; 0n 

W252. SEER 

1236 =Im= 

1240 INPUT "ENTER INPUT MODE (@=34974,1=FILE) O=DEFAULT”,Im 
1244! 

1248! INPUT MODE ZERO IS FROM THE DATA AQUISITION SYSTEM 

12525 1F Im=0) (HEN 


1256 BEEP 

1268 INPUT "GIVE A NAME FOR THE RAW DATA FILE" ,O2files 

1264 PRINT USING “16X,"“File mame: “",14A"iD2files 

1268! 

1272) CREATE BDAT FILE ON THE MASS STORAGE MEDIA 

1276 CREATE 8DAT D2file$,60 

1280! CREATE AN INPUT/OUTPUT LINK TO OPEN FILES 

1284 ASSIGN @File2 TO D2file$s 

1288! 

jes aa CREATE DUMMY FILE UNTIL Nrun KNOWN 

1296 DifileS="DUMMY" 

1302 CREATE 8DAT Dlfile$,60 

1304 ASSIGN @Filel TO Dlfiles 

1308 OUTPUT @FileliDates 

Pai! 

1316! CREATE A PLOT FILE 

1320 8EEP 

1324 INPUT "GIVE A NAME FOR THE PLOT FILE" ,Pfiles 

1328 CREATE BDAT Pfilet,30 

bSa2 ASSIGN ®Plot TO Pfiles 

1336 BEEP 

1340! 

1344! IOTC = NUM8ER (TOTAL) OF DEFECTIVE THERMOCOUPLES 

1348 INPUT “ENTER NUMBER OF DEFECTIVE TCS (@=*DEFAULT)*,Idtc 
WSS ! LOTC = LOCATION OF DEFECTIVE THERMOCOUPLE 

1356! 

1360 IF Idtc=@ THEN 

1364 PRINT USING “16X,""No defective TCs exist”"”* 

1368 ELSE 

13si2 PRINT USING “16X,""Defective Thermocouples Indicated by -99.99""" 
1376 END IF 

1380! 

1384 8EEP 

1388! DEFECTIVE THERMOCOUPLES MAY 8E IN CHANNELS 40-69 

1332! THEMOCOUPLES ARE ENTERED AS DEFECTIVE 8Y COMPUTER CHANNEL NO. 
1396! JOTC=COUNTER IN LOOP FOR DEFECTIVE THERMOCOUPLES 

14QQ! 

1404 IF Idtc*®@ THEN 

1408 FOR Jdtc=@ TO Idtc-1 

1412 INPUT "ENTER DEFECTIVE TC LOCATION (BY COMPUTER CHANNEL NUMBER ) 
Teoadtetidecs) 

1416 BEEP 

1420 NEXT Jdtc 

1424 END IF 

1428 PRINTER IS 701 

1432 OUTPUT @Filelitdtc(+#) 

1436 !, 

1440 ! Im=1 option (THIS OPTION ALLOWS DATA ENTRY WITH DATA FILE? 
1444 ELSE 

1448 8EEP 

1452 INPUT "GIVE THE NAME OF THE EXISTING DATA FILE” ,D2files 
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1456 
1460 
1464 
1468 
1472 
1476 
1480 
1484 
1488 
1492 
1496 
1580 
1504 
1508 
Lei 
1516 
1520 
1524 
Boze 
Po3e 
1546! 
154@ 
1544 
1548 
PSS2 
1556 
1562 
1564 
1568 
Sy gp 
1576! 
158 
1584 
1588 
1S32 
P5356! 
1620 
16Q4! 
16@8 
1652 
1616 
162@ 
1624 
1628! 
1632 
1636! 
1642 
1644 
1648 
1652 
1656 
166@! 
1664! 
1668 
1672 
1676 
1682 
1684 
1688! 
less 


PRINT USING “16X,""File name 
ASSISN ®File2 TO D2files 
ENTER 8File2:iNrun 

ENTER @File2;:Dold$,Ldtc(*),Itt,Bop,Nht ,Natp,Nrt ,Corr 
BEEP 

INPUT "GIVE A NAME FOR PLOT FILE" ,Pfile$ 

CREATE 8DAT Pfile$,30 

ASSIGN Plot TO Pfile$s 


"" 114A"; O2fi les 


PRINT USING "16X,""This data set taken on: °",14A",Dold$ 
BEEP 
PRINTER IS 1 
PRINT USING “4X,"““SELECT TU8E TYPE"”" 
PRINT USING "6X,""@ SMOOTH" "“ 
PRINT USING “6X,""“1l FINNED(1S/IN) """ 
PRINT USING “GX,”° 2 HIGH FLUX °*"* 
PRINT USING “6X,°"3 TUR8O-B “"" 
INPUT Itt 
END IF 


IF Im=1 THEN GOTO 1768 
PRINFER- IS 4 


IF Im=@ THEN 


PRINT USING "“4X,""Select tube type"”” 
PRINT USING “6X,°" @ Smooth “"“ 
PRINT USING “6X,"" 1 FINNED 19/IN (DEFAULT)"*" 
PRINT USING "6X,"" 2 HIGH FLUX"°*" 
PRINT USING "6X,"" 3 TURBO-B""" 
PRINT USING “6X,"" 4 GROWTH" "" 
PRINT USING "6X,"" S GROWTH""® 
PRINT USING "6X,"" 6 GROWTH""" 
ITT=TUBE TYPE 
INPUT Itt 
OUTPUT @Filel: Itt 
END IF 
PRINTER IS 7@1 
Itt=2 
PRINT USING "16X,""Tube Type: eee OM ey lett 
BEEP 
8op=@ 


INPUT "ENTER 8ULK OIL % 
OUTPUT @Filel:8op 
PRINT USING “16X,° "Bulk o1l’e” 0D" iB8op 


(DEFAULT#=@%) “,Bop 


BEEP 

NHT=NUM8ER OF HEATED TUBES 

Nht=5 

INPUT “Enter number of heated instrumented tubes(default=5)" ,Nht 


OUTPUT @FileliNht 
PRINT USING “16x,° 
SEEC 


Natp=Number of active dummy pairs 

Natp=2 

INPUT “Enter number of active dummy pairs (Default=@)",Natp 
OUTPUT ®8Filel;Natp 

PRINT USING “16X,""Number of active dummy pairs="",DD".Natp 
8EEP 


NRT=NUM8ER OF ADDED HEATED TU8ES TO ENHANCE 8UNDLE EFFECT 
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“Number of heated instrumented tubes="",D0D":Nht 


1696 Nrt=2 

17@0 INPUT "Enter number of added heated tubes from simulation heaters(Default: 
Q@)",Nrt 

1704 OUTPUT ®@FileliNrt 

1708 PRINT USING “16X,""“Number of added heated tubes( from simulation heaters )=“ 
*.O0D7. Net 

17iZ2 BEEP 

1716! 

172@! CORR IS CORRECTION FOR INSTRUMENTED TUBE HEIGHT 

1724 Corr=0 

1728 INPUT “WANT TO CORRECT TSAT FOR TUBE HEIGHT (@=YES(DEFAULT),1=NO)",Corr 


1752 IF Corr=@ THEN PRINT USING “16X,""TSAT 18 corrected instrumented heat 
ed tube height°"" 
1736 IF Corr=1 THEN PRINT USING "16X,""“TSAT 185 NOT corrected for instrumen 


ted heated tube height""” 

174Q@ OUTPUT @FileliCorr 

1744 8EEP 

1748! ILQU*INPUT MODE: LIQUID, VAPOR,OR LIQUID VAPOR AVERAGE 
1752 Ilav=2@ 

1756 INPUT "SELECT (@=LIQ( default), 1=VAP,2=(LIQ+VAP)/2)" ,Ilav 
176@! 

1764! DlA=Diameter at thermocouple positions (meters) 
1768 DATA .@122,0.0098,0.0106,0.0116,0,0,@ 

1772 READ Dla(+#) 

1776 Dl=Dla(Itt? 

178Q! 

1784! D2=Diameter to base of fins (outside dia for smooth)(meters ) 
1788 DATA .@158,8.0125,0.0158,0.01415,20,0,0 

1792 READ D2a(*#) 

1796 D2=DZa( Itt) 

18@Q! 

1804! Di=Inside diameter of unenhanced ends (meters) 
1808 DATA .@132,0.0109,0.0116,0.0127,0,0,2 

1812 READ Dia(t»#) 

1816 DiszDia(Itt) 

182Q! 

1824! Do=OQutside diameter of unenhanced ends (meters) 
1828 DATA .@15875,0.0125,0.015875,0.01415,0,0,2 

1832 READ Doa(*) 

1836 Do=Doa(Itt) 

184Q! 

1844! L=Length of enhanced surface (meters) 

848 DATA® 22032 ,.2032,.20352,.2032,.2032,.2032,.2032 
18S2 READ La(#) 

1856 Le=La(Itt) 

186Q! 

1864! Lu=CORRECTED Length of unenhanced surface at the ends (METERS) 
1868!LU=LFIN + THICKNESS/2 

1872 DATA .@261, .0254,.0264,0.0258,0,0,0 

1876 READ Lual+) 

188@ Lu=Lua( Itt) 

1884! 

1888! LV=corrected length of 3 inch finned like end 
1892 DIM Lva‘6) 

1896 DATA .@769,.0762,.0773,0.0766,0,0,0 

19@0 READ Lva(+*) 

1904 Lv=Lva( Itt) 

1908! Kcua=Thermal Conductivity of tube 

1912! DATA 401,20,2@,%,9,0,2 

1916! READ Kcua(#) 
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1922! Keu=kKcual( ltt) 

1924 A=PI#(Do0%2-D1°2)/4 

1928 P=PI+«Do 

1932 J=i 

1936 Sx=@ 

194@ Sy=@ 

1944 Sxs=@ 

1948 Sxy=@ 

1952 Repeat: ! 

1956! 

1960 IF Im=@ THEN 

1964! Otld=desired temperature of liquid 

1968 Dtld=47.S UR= 113 

1969! Dtld=2.2 IR-114 

a7 2 Ido=2 

1376 ON KEY @,15 RECOVER 1952 

198@ PRINTER IS 1 

1984 PRINT USING “4X,""SELECT OPTION “"" 

1988 PRINT USING “6X,"““@=TAKE DATA“"" 

ES is Fd PRINT USING “6X,""1=SET HEAT FLUX""”" 

1996 PRINT USING “6X,""“2=SET Tsat (DEFAULT SET FOR R-113)""" 

2282 PRINT USING "4X,""NOTE: KEY @ = ESCAPE""" 

2004 ! Ido=desired option 

2008 BEEP 

2Q12 INPUT Ido 

2016! 

2020 BEEP 

2024! Set default value for input 

2028 IF Ido>2 THEN Ido=2 

Z0352! Take data option 

2036 IF Ido=@ THEN 2442 

2240! 

2044! LOOP TO SET HEAT FLUX (FOR TOP INSTRUMENTED TUBE) 

2048 IF Ido=1 THEN 

29S2 Dadp=100020 

2@S6 PRINT USING “4X,""Qdp QDPsim Nrt Qdpaux 
Q@tot""* 

2068 PRINT USING “4X,""“(W/m%2) (W/m*2) (W/m*2) 
cWpe” 

2064 Err=] 

2@68 ! Reset,read channel 25-3@,automatic scaling 

2072! Channel 25*au~ amps,26=sim amps,27#inst volts,28=sim volts,29=aux 

volts,3@=i1nst amps 

2076 OUTPUT 7@9; "AR AF2S ALS@ VRS" 

208@ FOR I=1@ TO 11 

2084 OUTPUT 709; "AS SA" 

2088 ENTER 709: Amp( 1) 

2092 NEXT I 

2096 FOR I=@ TO 2 

2100 OUTPUT 709: "AS SA" 

2104 ENTER 709;Volt(I) 

2108 NEXT I 

a ns OUTPUT 709: "AS SA" 

AIRS ENTER 709; Amp(Q) 

212Q! Calculate actual heat flu» 

2124 , Q(@)=60#Volt(O)#*Amp(@) 

2r2o Qdp(@)=Q(@)/(PI*O2eL ) 

7s le Qsim=62#20@*#Volt(l d#Amp( 11) 

ualeis Qdpsim=Qsim/(PI*#D2*.2032+5 ) 

2142 Qaux»=60*20@*#Volt( 2 #Amp( 1) 
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2144 Qdpaux=Qaux/(PI*.Q@160«.177&8«4) 


2148 Qtot=Q(0)«Nht+QsimtQaux 
2USZ Nrt=Qdpsim/Qdp(Q0) 
2156 IF ABS(Aqdp-Dadp)>Err THEN 
2162 IF Aqdp>Dadp THEN 
2164 BEEP 4000,.2 
2168 EESE 
PAM I BEEP 3250.22 
2176 END IF 
2180 IF Nrt<.1 THEN Nrt=@ 
2184 IF Qdpaux<10@ THEN Qdpaux=0 
2188 IF Qdpsim<100 THEN Qdpsim=Q 
2132 PRINT USING “4X,2(MZ.3DE,2X),2X,(MDD.DD),2X,2(NZ.3DE,2Xx ) ada 
(Q@),Qdpsim,Nrt,Qdpaux,Qtot 
2196 WAIT 2 
2200 GOTO 2076 
2204 END IF 
2208 END IF 
2212! 
2216! LOOP TO SET Tsat 
2220 IF Ido=2 THEN 
2224 IF Ikdt#=1 THEN 2240 
Z22ec BEEP 
22oc INPUT "ENTER DESIRED Tsat (DEFAULT=#=47.5 C - R-113)",Dtld 
2255" INPUT “ENTER DESIRED Tsat (DEFAULT#=2.2 C - R-114)",Dtld 
2256 Ikdt=1 
242 Oldl=0 
2244 O1ld2=0 
2248 Nn#=] 
2252 Nes=#=Nn MOD 15 
2256 Nn#=Nn+1 
2260 IF Nrs#l1 THEN 
2264 PRINT USING "4X,"" DTsat Tldl Tldz Tlbb 
Tvat Tvab Tlav ais 
2268 END IE 
222! Read thermocouple voltages for vapor, liquid 
2276 OUTPUT 709; "AR AF@ ALS VRS" 
2280! Sample each thermocouple 20 times and report temp for each the 
rmocouple, vapor#0,1,2; liquid=384 
2284 FOR I=@ TO S 
2288 Sum=Q 
2292 OUTPUT 709; "AS SA” 
2296 FOR Ji=l1 TO 20 
2380 ENTER 709;Eliq 
2304 Sum=Sumt+E1lig 
2308 NEXT Ji 
PES 4 Emf( I )=Sum/20 
eol6 TCI =FNTvsv(Emf(1)) 
2o2 NEXT I 
eac4 | Compute average temperature of liquid 
fis Gia: Tlav=(T(3)+T(4))#.5 
eoce' Compute average temperature of vapor 
2335 Tvavl=(T(@)4T(1))/2 
2535 Tvave=T(2) 
2556 Tvav=(T(Q0)+T(1)4T(2))73 
2348 : IF ABS(Tlav-Dtld)>.2 THEN 
2344 IF Tlav>oDtld THEN 
2348 BEEP 4000,.2 
Zoo ELSE 
2356 BEEP 250,.2 
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2562 ENO IF 


2364 ELSE 

2368 IF ABS(Tlav-Otld)>.1 THEN 
Zale IF Atld>Otld THEN 
es76 BEEP 3000,.2 
2380 ELSE 

2384 BEEP 822,.2 

2388 ENO IF 

2592 ENO IF 

2396 ENO IF 

2402 Errl=Tlev-Old! 

2404 Oldl=Tlav 

2428 Err2=Tvav-Old2 

2412 Old2=Tvav 

2416 PRINT USING “4X,7(M000.00,3X)";Otld,7T(3),7T(4),T(S), Tvavl, Tvav2,Tlea 
v 

2422 WAIT 2 

2424 GOTO 2252 

2428 END IF 

2432! 

2436! TAKE DATA IF Im=@ LOOP 

2442 IF Ikol=1 THEN 2452 

2444 BEEP 

2448 Ikol=1 

2452 OUTPUT 789; “AR AF@ ALS VRS" 
2456 FOR I=@ TO S 

2462 OUTPUT 709; “AS SA” 

2464 Sum=2 

2468 FOR Ji=1 TO 2@ 

Fal ee ENTER 7Q9;E 

2476 Sum=Sum+& 

2482 IF I>2 THEN Et(Ji-1l)=#€ 
2484 NEXT Ja 

2488 Kd1=Q 

2492 IF I>2 THEN 

2496 Eave=Sum/20 

2500 Sum=@. 

2504 FOR Jk=@ TO 19 

2508 IF ABSC(Et( Jk )-Eave)<S.@E-6 THEN 
252 Sum*™SumtEt (Jk) 
2516 ELSE 

2520 Kdl=Kdl+1 
2524 ENO IF 

2526 NEXT Jk 

fils fe fe IF I>2 THEN PRINT USING “4X,""Kdl = “*,OD";Kdl 
2536 IF Kdl>1@ THEN 
2542 BEEF 

2544 BEER 

2548 PRINT USING “4X,""Too much scattering in data - re 
peat data set""" 

age 650TO 1982 

2556 ENO IF 

2562 ENO IF 

2564 Emf( I )=Sum/( 2@-Kdl ? 

_568 NEXT I 

Ere OUTPUT 709; "AR AF4@ ALES VRS" 
2a76 FOR I=6 TO 3S 

2582 OUTPUT 709: "AS SA” 

2584 Sum=@ 

2958 FOR Jiz=l TO S 
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Soe 
2546 
2620 
-604 
2628 
eos! 
2616! 
2622 
2624 
2640 
26S. 
2636 
264Q! 
2644! 
2648 
26s2 
2656 
2662 
2664 
2668! 
2672 
2676 
2680 
2684 
2688 
2692 
2696 
2700 
27Q4 | 
2798 | 
2tu2 
27165 
ened 
2724 
2728 
Zoe 
2756 
2742 
2744! 
2748 
Pah om Fe 
2290 
-76Q! 
2764 
Cree 
Ra 7 
ceio! 


ENTER 7Q9;E 
Sum=Sum+& 
NEXT Ji 
Emf( 1 )=Sum/S 
NEXT I 


READ VOLTAGES (27=Inst ,28=Sim,29=Aux ) 
OUTPUT 70S; GR AFZ7 ALZS VRS- 
FOR I[=8@ TO 2 
CUTPUI 703; As sa. 
ENTER 709;Volt( 1) 
NEXT I 


READ CURRENTS (30-34=Inst tubes: 3S-39=ACTIVE Dummy ) 
OUTPUT 7209; "AR AF32@ AL39 VRS" 
FOR I=@ TO 9 
GUTPUT 709:°AS SA” 
ENTER 709; AmpC I) 
NEXT I 
Read Currents(2S=Aux amps,26=Sim amps) 
OUTPUT 709; "AR AFZ2S ALZE VRS" 
FOR I=1@ TO 11 
OUTPUT 7093 As SA 
ENTER 7209; Amp( 1) 
NEXT I 
EESE 
ENTER @File2;sEmfl*) Volt (*#),Amp(*) 
END IF 


CONVERT EMF’S TO TEMP, VOLT,CURRENT 
FOR I=Q2 TO 35 
TCI )=FNTvsv(Emf(.I)) 
IF I>S AND Idtc>Q@ THEN 
FOR Ii=Q@ TO Idte-l 
IF Udtc( Ii =I -S+39 THEN T(I )#-99.99 


NEXT Ti 
END IF 
NEXT I 
Ntc=nr of thermocouples 
Ntc=6 


FOR Tix@ TO 4 
Q(T1 )=6Q*Volt(Q)*Amp( i) 
Twa=Average temperature of the wall 
Twa( Ti )=Q 
Ndtc=2 
FOR T=1 TO Nte 
Nn is counter in temp array, start at 5 ‘this 1s the 


ermocouple in the tube ban ) 


2780 
2784 
2788 
oese 
2796 
~ 822 
2804 
2808 
2812 
-816 
2820 
a2 ea | 


om we oe 


Nn=I1#64I+5 
IF ABS(T(Nn)) 99 THEN 
T(Nn > =-99.99 
Ndtc=Ndtctl 
ELSE 
Twa( Ii )=Twa( Ii »+T(Nn) 
END IF 
NEXT I 
Twa( ti d=Twa(T1)/(6-Ndtc ? 
NEXT I1 
Tlav=(T(3)+T(4))/72 
Tvav=T(2) 
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first th 


2824! 
aa | 
Zac! 
2852 
2950 
2840! 
2844! 
2848 
Fl 8 bo 
2956 
2860 
2864 
2868 
20720 
2876 
2880 
2884! 
2888! 
2892 
2896 
2980 
2904 
2908 
gole 
2916 
2420 
2924 
2923 
ego! 
~936! 
2948 
2944 
2948 
Pid jo 
2396 
2960 
2964 
2968 
E02 
2976 
2980 
2984 
2988! 
2992 ! 
2336 
3000 
3004 
3008 
S@1c 
3016 
3020 
3024 
3028 
3032 
5036 
5042 
5044 
3048 
3052 
3056 


Fvav=( FOO)+1 (1 i+Ie prs 


Tlav=T( 


>) 


Tcu=Twa(Q@) 


Kcu=FNK 


Twf IT )=Twall )-Q¢ 1 )*LOG(D2/D1)/(2*PI*Kcust ) 


iE 
Le 
Us 
IF 
hE 
TF 
NEXT I 


COMPUTE 


Tfilm=(Tw(@)+Texs)#.5 


Rho=FNR 
Mu=FNMu 
K=FNK(T 
Cp=FNCp 
8eta=FN 
N2=Mu/R 
Alpha=K 


eut leu) 


ITlqgv=@ THEN 
Tlqv=1 THEN 
Ilqv=2 THEN 
Corr=1 THEN 
Corr=@ THEN 
Corr=@ THEN 


THERMAL CONDUCTIVITY OF GOPFER 


MTP s CURVE FIT NOT AUAIE USE ARRAY KCUX*) 
FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
FOR I=@ TO 4 


Texs=Tlav 
Texs=Tvav 


Texs=(Tlavt+T(2))+* 
Thetab( I )=Tw( I )-Texs 

Thetab( I )=Tw( I )-( Texs+.Q@56+I1*.129) 
Thetab( I )=Tw( I )-( Texs+. 054+] #.144) 


VARIOUS PROPERTIES 


ho(Tfilm) 
(Tfilm) 
film) 
CTfaim) 
Seta(Tfilm) 
ho 
/(Rho#Cp ) 


Pr=Ni/Alpha 


COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 


FOR UNENHANCED END(S) 


Lu=Lua( 


Itt? 


Hbar=190 
Fe=(Hbar+*P/(Kcut#A))*.Se*Lu 


Tanh=FN 


Tanh( Fe) 


Theta(@)=Thetab(@)*Tanh/Fe 
Xx~=(9.81*Beta#Thetab(@)*D0°3*Tanh/(Fe#Ni*Alpha))*.166667 


Yy=(1+(.559/Pr )*(9/16))*( 8/27) 


Hbarc=K/Do#( .6+.387#Xx/Yy )*2 


IF A8S( (Hbar-Hbarc )/Hbarc )>.@@01 THEN 


Hbar=(HbartHbarc)#.5 


60T 
END IF 


COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S ) 
Q1(@)=( Thetab(@)*#Tanh)«#( (Hbar*P#kKcutA)” .5) 


Qq=Q1(0 
Z=Z+1 
IF Z=1 


0 2948 


1+0Q 


THEN 


Lu=Lv 


GOTO 
ENO “fF 
Z=Q 


~944 


Qipct=Qq/0(0) 


Qq=0 


As=PI*#D24L 


FOR I1= 


@ TO 4 


Olt trI=Oiect*Qt 12) 


Odpet li o=(GC11)=G1(11007As 
Htube(I1 )=Qdp(I1)/Thetabil.) 


NEXT [1 
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A 


tELEM TEMPERATURE 
IDENSITY 

I'VISCOSITY 

' THERMAL CONDUCTIVITY 
VSPEGIE LC cneAT 
ITHERMAL EXPANSION 
IRINEMATIC VISCOSITY 
THERMAL DIFFUSIVITY 
'PRANDTL 


IR-114 
IR-113 


3062 
5064! 
3@68 ! 
3@72 
3076 
3280 
3284 
3088! 
5@92! 
3096 
3180 
3104 
31@8 
Eo UF 
Tvav 
3116 
3122 
3124 


Shes 


PRINTER [S 7@1] 


RECORD TIME OF DATA TAKING 
IF Im=@ THEN 

OUTPUT 7809; “TD” 

ENTER 709: Told$ 
END IF 


OUTPUT DATA TO PRINTER 

PRINTER Ts e701 

PRINT 

PRINT USING “1@X,""Data Set Number = "",00D,2X,14A";J,Tolds 

PRINT 

PRINT USING "10X,"" Tvl Tv2 Tv3 {alll itez Its 


Tldav """ 


PRINT USING "10x ,8(MDD.DD, 2x d°; T(0),TC1),702),7(3),104),1(S), iva 
PRINT 


PRINT USING "6X,""“Tube Wall Temperatures (Deg C)? Tnave Qdp 
Thetab"“"" 
PRINT USING "6X,""8 l 2 3 4 5 6 (Deg C) (W/m* 


CIO MeN Zeke) i) ae 


Blog 
3136 
3142 
3144 
3148 
S1s2 
5156 
316@ 
3164 
3168 


3172 


J 3=@ 
FOR 11=@ TO Nht-1l 
FOR Ji=@ TO S 
To(J1)=T(11*5+Jj+6) 
Jy=Jyt+l 
NEXT Ja 
Jjy=litl 
FOR Ji=@ TO 4 
resp lt). 
NEXT Ja 
PRINT USING “6X,0,1X,7(MD0.00),1X ,2¢M%Z.3DE),1X,1(MDD.00) "; Tne 


o(@),Tp(1),7Tp(2),Tp(3),7Tp(4),To(S), Twal li), Qdp(1l1),Htube(I1),Thetab(Ii) 


3176 
5180 
3184 
3188 
o1g2 


tors fa 


NEXT Ia 
Ok=1 
IF Im=@ THEN 
BEEP 
INPUT “OK TO STORE THIS BATA SET (1#Y( default ),@=N)7", 0k 
END IF 
J=the counter for data sets 
IF Ok=1 OR Im=1 THEN J=J+1 
IF Ok=1 AND Im=@ THEN OUTPUT @FileliEmf(#),Volt(#),Amp(*#) 
IF Im=1 OR Ok=lL THEN OUTPUT ®@Plot:iQdp(*),Htube(*), Thetab( *) 
Go_on=1 
IF Im=@ THEN 
BEEP 
INPUT "WILL THERE BE ANOTHER RUN (1=Y‘ dafault),@=N)?",Go_on 
Nrun=J 
IF Gc_on=@ THEN 3272 
IF Go_on-.*@ THEN Repeat 
EUSE 
IF J Nrunt! THEN Repeat 
ENO IF 
St=1 
BEEP 
INPUT "ARE YOU SURE YOUR READY TO TERMINATE (1=Y(DEFAULT),@=NQO)?7" St 
Go_on=1 
IF St°@ THEN 35280 
IF St=@ THEN GOTO 3242 
IF Im=@ THEN 
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3284 BEEP 


3288 PRINT 

Goede PRINT USING “10X," “NOTE: ~",Z2Z.,°" d@ta cuns weré stored in file 
, 190A"; J-1,D02files 

3296 ASSIGN @Filel TO * 

3320 OUTPUT @File2;Nrun-1 

3304 ASSIGN ®Filel TO Difiles 

3308 ENTER @Filel;:DateS,Ldtc(#),Itt,Bop,Nht ,Natp,Nrt,Corr 

Bole OUTPUT @FileZ2:Date$,Ldtc(+),Itt ,Bop,Nht ,Natp,Nrt,Corr 

S516 FOR I=! TO Nrun-l 

3320 ENTER ®@FileliEmf(#), Volt(#),Amp(*#) 

3324 OUTPUT @File2:Emf(#),Volt(*),Amp(«#) 

S42 NEXT I 

S552 ASSIGN @Filel TO * 

3336 PURGE “DUMMY " 

3340 END IF 

3344 BEER 

3348 PRINT 

Z552 PRINT USING "10X,""NOTE: °",ZZ,"" X~Y pairs were stored in plot data f 
tle "" 10A"; J=1,Pfiles 

3356 ASSIGN @File2 TO * 

3362 ASSIGN @Plot TO * 

3364 BEEP 

3368 SUBEND 

aoc 

3376! XXXXXXXXXXXXXXXXXXXXKXKXKXXKXKXKXAKXKXKKKKAKXKKXKXKKKKXKKK XK KKK KKK KKK KKK KKK K KKK KKK KKK KKK 
338! 


3384 DEF FNKceu(Tcu) 
3388! OFHC COPPER 


5392 Tk#Teut273.1S iG TOK 

3396! Keu=434-.112+4Tk 1250-3@0K USE FOR R-114 62.2 C 
3400 Keu=433.0-.14#Tk '200-4@0K USE FOR R-113 047.5 C 
3404 RETURN Kceu 

3408 FNEND 

6412! 


5416 DEF FNMu(T) 
S420! CURVE FIT OF VISCOSITY 


3424 Tk=T+273.15 'C TO K 

3428! Mu=EXP(-4.4636+(1011.47/Tk))#1.@E-3 '!R-114 170-3560 K 
3432 Mu=.0000134#(10*(SO03/(Tk-2.15))) IRid3 

3436 RETURN Mu 

344@ FNEND 

3444! 


3448 DEF FNCp(T) 

3452! CURVE FIT OF Cp 

3456 Tk=T#273.15 Le mOeK 

3460! Co=.40188+1 .GS@Q@7E-S#Tk+1.51494E-6eTh *2-6.678S53E-10+Tk*°3 IR-114 180-400 k 
3464 Cp=(929+1.03+*T)*.0@1 IR-1135 

3468 Cp=Cp+1Q02 

5472 RETURN Cp 


3476 FNEND 

5480! 

5484 DEF FNRho!T) 

9463 (TReTt2/3.15 'C TO K 


S492 Meml-(1.0*Tt 7753.95) 'K TO PR 

9496!) Ro=-36.52461.146414 9%" 1/3)4+16.41801SeX+17.476858ex .S+1.1198l0e*x*2 
3500! Ro=Ro/ .062428 b= 1 14 

G04) ho-l75207479E+5-14(2. 0 1S6546+7*. 35782 91E-35) a loaae 22S. 

3508 RETURN Ro : 

S512. FNEND 
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SSL50 


5.0 
3524 
5528 
asoe 
35956 
3540 
3544 
3548 
g5o5 
3996 
556@ 
3564 
S560 
Sore 
S575 
358 
3584 
3588 
Sone 
2596 
5608 
3604 
3608 
able 
3616 
362@ 


3624! 


3628 
S632 
3636 
3640 
3644 
3648 
Sbo2 
3656 
3662 
3664 
3668 
3672 
3676 
5680 
3684 
3688 
6692 
2526 
3702 
3704 
3708 
S7)2 
S715 
orcO 
3724 
oe 5 
Beye 
5736 
3742 
5744 
3748 
Se 


DEF FNPriT) 1'GOQOD FOR R-114/R-113 
Pr=FNCp(T)*FNMu( T)/FNK(T) 

RETURN Pr 

FNEND 

DEF FNK(T) 


T< 360 K WITH T IN C 
K=.071-.@@@261+#T 
RETURN K 

FNEND 


DEF FNTanh( Fe} 
P=EXP( Fe) 

Q=EXP( -Fe) 
Tanh=(P-Q)/(P+Q) 
RETURN Tanh 
FNEND 


DEF FNTvsvi(V) 
COM 7Ge7ti7)) 
T=C(Q) 

FOR I=l1 TO 7 
T=T+C( I )#U*] 
NEXT I 

RETURN T 
FNEND 


DEF FNBeta(T ) 
Rop=FNRho( T+.1) 
Rom=FNRho(T-.1 ) 
Beta=-2/(RoptRom)#*(Rop-Rom)/.2 
RETURN Beta 
FNEND 
DEF FNPoly(X) 
COM /Cply/ AC1@,10),C0(10),8(4),Nop,Iprnt,Opo,llog 
X 1=X 
Poly=B(@) 
FOR I=1 TO Nop 
IF Tlog=1 THEN X1=LOG(X) 
Poly=Poly+B(I )#X1°I 
NEXT I 
IF Tlog=1 THEN Poly=EXP( Poly) 
RETURN Poly 
FNEND 


EXXXXXMXXXXXXXXXAXXXKXKXKXAXAXX KKK KKK K KKK KK KKK KK XXX KK XXKXKX KK XXX KXKXXKXKXKKKKXKXAKXXKXKK 


SUB FPoly( Of ples C*) Np. lin) 

DIM ORi105, S010). Sy0 120554 129 xs 1000, Yy0 100). vel 
COMPVEply? At 1OSI1G °C 105.604) No iprnt.0po, log 
COM Vs kx yy? pts). Ypc 

FOR 1=@ TO 4 

B( 1 )=Q 

NEXT I 

Im=1 

BEEP 

INPUT “ENTER DATA FILE NAME" ,Ofrles(Q) 

BEEE 

INPUT “ENTER NUMBER OF X-Y PAIRS" ,Np 
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BEEP 
INPUT "LIKE TO EXCLUDE DATA PAIRS (1=Y,Q=N( DEFAULT: )?",Ted 
IF ITed=1 THEN 

BEEP 

INPUT "ENTER NUMBER OF PAIRS TO 8E EXCLUDED", Ipe~ 
END IF 
ASSIGN @File TO DfirleS(Q) 
N=2 
BEEP 
INPUT "ENTER THE ORDER OF POLYNOMIAL (DEFAULT=2) ",N 
FOR [=@ TON 

Sy(1I)=@ 

S«( 1 )=@ 

NEXT I 
IF ITed=1 AND Im=1 THEN 

FOR I=1 TO Ipex 

ENTER @File:Xy(*#) 

NEXT I 
END IF 
FOR I=1 TO Np-Ipe~ 

ENTER @FileiXy(*#) 

IF Opo=@ THEN 
Y=Xy(Itn-1) 
X=Xy(lititn) 

END IF ; 

IF Opo=1 THEN 
Y=Xy(StItn) 
X=Xy(11+Itn) 

END 1F 

IF Opo=2 THEN 
Y=Xy(S+Itn) 
X=Xy(Itn-l) 

END IF 

IF [Tlog=1 THEN 
X=LOG(X) 

Y#LOG(Y) 

END IF 

Xx( J )=X 

Yy(I)=Y 

R(Q)=Y 

Sy(@d=Sy(QI+Y 

S¢(1)=X 

S*(1=Sx0 1 )4+X 

FOR J=1 TON 

R( J )=R(J-1 )#X 
Sy(Jd=Sy( J +RC I) 
NEXT J 
FOR J=2 TO Ne2 
SCJ =S( J-1 #X 
Sa6 J 0=SxO J 450d) 
NEXT J 
NEXT I 
S«(Q)=Np 
FOR I[=@ TON 
C(I )=eSyf 1) 
FOR J=@ TON 
AC I,J )=Ss.(I+J) 
NEXT J 
NEXT I 
FOR I[=@ TO N-1 
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3996 CALL Divide(l) 
4000 CALL Subtract(I+1) 
4004 NEXT I 

4008 B(N)=C(N)/A(N,N) 

4012 FOR I=0 TO N-1 


4016 BCN=1=1 }=C(N=1-T > 

4220 FOR J=@ TO I 

4024 BON=1=))=8¢N=)=f J-ACN-1=1 (N=) eBCN-2 > 
4@28 NEXT J 

4032 BC N-l-Iv=BCN=1=1/ACNS TSP N=1=1) 


4036 NEXT I 
4040 IF Iprnt=@ THEN 


4044 PRINT USING “12X,""“EXPONENT COEFFICIENT""” 

4048 FOR I=@ TON 

4052 PRINT USING “15X,00,5X,MO.70E”";1I,8¢(1) 

4056 NEXT I 

4Q62 PRINT " “ 

4064 PRINT USING "12X,""DATA POINT X Y Y(CALCULATED) OI 
SCREPANCY""“ 

4268 FOR I=1 TO Nop 

4072 Yo=B(@Q@) 

4076 FOR J=1 TON 

4Q@8Q YerYo+B( J )#Xx(I)°J 

4084 NEXT J 

4088 D=Yy(I)-Ye 

4092 PRINT USING “15xX,30,4X,4(MD.S50E,1X) vex Gly) Youn 
4096 NEXT I 

4100 END IF 


41@4 ASSIGN @File TO # 

41@8 SUBEND 

4112! 

4116 SUB Divide(M) 

4120 COM /Cply7 A 10,10),C(10),B8(4).,N, locnt ,Op0, 1 los 
4124 FOR I=M TON 


4128 Ao=A(I,M) 

4132 FOR J=-M TO N 

4136 ACI, J )#A(I,J)/Ao 
4140 NEXT J 

4144 C(I )=C(1)/Ao 


4148 NEXT I 
4152 SUBEND 
4156 
4160 SUB Subtract(K) 

4164 COM /Coly/ AC(1@,10),C(1@),B8(4),N,Iprnt,Opo,liog 
4168 FOR I=K TO N 


— 


4172 FOR J=K-1 TO N 

4176 AOS ACK ly sm ola? 
4180 NEXT J 

4184 Cty= Chl -eCL) 


4188 NEXT I 

4192 SUBENO 

4196! 

42@0 SUB Plin 

4204 COM /Cply/ AC1@,1@),C(1@).8(4),N,Iprnt,Opo,llog 
4208 COM /Xxyy/ Xx«(S),Yy(S) 

4212 PRINTER IS 70S 

4216 BEEP 

422Q@ INPUT “SELECT (@=h/h@% same tube, l=h( HF )/h( sm)", Irt 
4224 BEEP 

4228 INPUT "WHICH Tsat (1=6.7,0=-2.2)",Isat 
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420 


doe 
4236 
4242 
4744 
4248 
e257 
4256 
4262 
4264 
4268 
Reve 
4276 
4280 
4294 
4288 
A292 
4296 
4300 
4304 
43208 
4312 
4316 
4320 
4324 
4329 
ASS. 
4336 
4340 
4344 
4348 
4352 
4356 
4360 
4364 
4368 
A372 
4376 
4380 
4384 
4388 
4392 
4396 
442Q 
4404 
4408 
4412 
4416 
44°22 
Aan 
4428 
ano 
4436 
4440 
4444 
4448 
4452 
4456 
4462 
4464 
4468 


Kmin=Q 

Kma«~=1Q 

Xstep=2 

IF Irt=@ THEN 

Ymin=Q 
Ymax=1.4 
Ystep=.2 
BES 
Ymin=0 
Ymar=15 
Ystep=5 

END IF 

BEEP 

PRINT “IN: SP1;IP 2300,2200,8300,5800; " 

PRINT “SC @,100,0,100:TL 2,Q:" 

Sfx=100/(Xmax-Xmin) 

Sf y=100/(Ymax-Ymin) 

PRINT "PU @,@ PD” 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#Sfx 
PRINT "PA";:X,",@:; XT; " 

NEXT Xa 

PRINT “FA 100.0; PU: ” 

PRINT "PU PA @,@ PD” 

FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
PRINT “PAO. 3Y.° YI” 

NEXT Ya 

PRINT "PA @,10@ TL @ 2” 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#Sfx 
PRENT “PAY SX 3 1005; XT” 

NEXT Xa 

PRINT "PA 100,100 PU PA 100,20 PD" 

FOR Ya=Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin)#Sfy 
PROENT "PO PA eIOO. ny 4 YE” 

NEXT Ya 

PRINT "PA 100,100 PU” 

PRINT. PA @.-2 Sk 1.5.27 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin)#Sfx 
PRINT “PA uk. 03 
PeorNin (esc. late ke, 

NEXT Xa 

PRINT “PU PA 0,0” 

FOR Ya=Ymin TO Ymax STEP Ystep 
IF ABS(Ya)<1.E-S THEN Ya=0 
Y=(Ya-Ymind#Sfy 
PRINT “PASO. #1.” 

PRINT ~ (OP <=4.--2 25:06 1Yau. 

NEXT Ya 

Xlabel$="011 Percent” 

[TP let=0 THEN 

YlabelS="h/hO%" 
» ESE 
YlabelS="h/hsmooth" 

ENO IF 

PRN cee eS. eo PA S010 EP 3s -LEN Alabel B72. 0326 

PRINT “PA -11,5@ CP Q@,";-LEN( YlabelS)/2+S5/6: "DI @,1:LB 
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7 AC eane ae 
omit batels.- ~ 


4472 PRINT “CP Q,2”" 

4476 Ipn=2 

448Q@ BEEP 

4484 INPUT “WANT TO PLOT DATA FROM A FILE (1=Y,0=N)?" ,Okp 
4488 Icn=Q 

4492 IF Okp#l1 THEN 

4496 BEEP 

4500 INPUT "ENTER THE NAME OF THE OATA FILE” ,0_fi les 
4504 BEEP 

4508 INPUT “SELECT (@=LINEAR, 1=LOG(X,Y)",Ilog 
4512 ASSIGN @File TO D_file$s 

4516 BEEP 

452@ INPUT “ENTER THE SEGINNING RUN NUMBER" ,Md 
4524 BEEP 

4528 INPUT “ENTER THE NUMS8ER OF X-Y PAIRS STORED" ,Npairs 
4532 BEEP 

4536 INPUT “ENTER DESIRED HEAT FLUX" ,Q 

4542 BEEP 

4544 PRINTER IS 1 

4548 PRINT USING “"4X,""“Select a symbol:"“" 

4552 PRINT USING “4X,""1l Star 2 Plus sipm. | 
45SG6 PRINT USING "4X,""3 Circle 4 Square""" 
4560 PRINT USING "4X,"°"S Rombus“"" 

4564 PRINT USING "4X,""6 Right-side-up triangle""" 
4568 PRINT USING "4X,°"7 Up-side-down triangle""" 
4572 INPUT Sym 

4576 PRINTER IS 7@S 

4580 PRINT “PU DI” 

4584 IF Sym=l THEN PRINT "SMe#" 

4588 IF Sym=2 THEN PRINT “SM+" 

4592 IF Sym=3 THEN PRINT "SMo" 

4596 Nn=4 

4600 IF Ilog=1 THEN Nn=1l 

4604 IF Md>1 THEN 

4608 FOR I=l1 TO (Md-1) 

4612 ENTER @File:Xa,Ya 

4616 NEXT I 

4620 END IF 

4624 Q12Q 

4628 IF Ilog=1 THEN Q=L0G6(Q) 

4632 FOR I=1 TO Npairs 

4636 ENTER @File:Xa,B(*#) 

4642 Ya=B(Q@) 

4644 FOR K=1 TO Nn 

4648 Ya=YatB(K )#Q*K 

4652 NEXT K 

46S6 IF ITlog=1 THEN Ya=EXP( Ya) 

466@ IF ITlog=@ THEN Ya=Q1/Ya 

4664 IF Irt=@ THEN 

4668 IF Xa=@ THEN 

4672 Yo=Ya 

4676 Ya=1 

4680 ELSE 

4684 Ya=Ya/Yo 

4688 ENO IF 

ASae ELSE 

4696 Hsm=FNHsmooth(Q,Xa,Isat ) 

4702 Ya=Ya/Hsm 

4704 END @IE 

4708 X«(I-1)=Xa 
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471. Yyti-l)=Ya 

4716 KX=(Xa-Xmin)#Sfy 

472@ Y=(Ya-Ymin)#Sfy 

4724 IF Sym>3 THEN PRINT “SM" 

4729 IF Sym<4 THEN PRINT “SR 1.4,2.4" 

A752, FRINGE “PA 2 xs Y, 

a7 500 1F Ssym23°. TRENSPRINT “SR -1.2,1.6~ 

£740 Sir Syms, THEN PRINT WCZ,4,99,0,-8,-4,8,0,824.,0.5— 
S744 mele  ovasooihen PRINT “WCS.0..39--3.-6.-3.6.5,0,5,-0° © 
6749 [6 Sym=o THEN PRINT “UCO,5.3,99..5,-8,76,0,3,8: ° 
ia. TE <Syme=7 THEN PRINT “UGO,-—5.3,99,-3,6.6,0.-5.-33 
4756 NEXT I 

476@ BEEP 

4764 ASSIGN &File TO * 

4768 ENO IF 

4772 PRINT “PU SM” 

A716 BEEP 

4780 INPUT "WANT TO PLOT A POLYNOMIAL (1=Y,@=N)?" ,Okp 
4784 IF Okp=1 THEN 


4788 BEEP 

4792 INPUT “SELECT (Q@=LINEAR,1=LOG(X,Y))",Ilog 
4796 Iprnt=l 

48Q@ CALLS PolyCitn) 

4804 FOR Xa=Xmin TO Xmamx STEP Xstep/25S 
488 Icen=Icntl 

4812 Ya=FNPoly( Xa) 

4816 Y=(Ya-Ymin)*#Sfy 

4822 X=(Xa-Xmin)#Sf~ 

4824 IF ¥<@ THEN Y=@Q 

4828 IF Y>1@@ THEN GOTO 4868 

4832 Pu=@ 

4836 IF Ipn=1 THEN Idf=Icn MOD 2 
4840 IF Ipn=2 THEN Idf=Icn MOD 4 

4844 IF Ipn=3 THEN Idf=Icn MOOD 8 
4848 IF Ipn=4 THEN Idf=Icn MOD 16 
48S2 IF Ipn=S THEN Idf=Icn MOD 32 
4856 IF Idf=1 THEN Pu=1 

4862 IF Pu=@ THEN PRINT "PA",X,Y,°PD" 
4864 IF Pu=1 THEN PRINT “PA",X,Y, "PU" 
4868 NEXT Xa 

4872 PRINT "PU" 

4876 Ipn=Ipnt+l 

4880 GOTO 448@ 

4884 END IF 

4888 BEEP 


4892 INPUT "WANT TO QUIT (1=Y,@=N)?",Iquit 
4896 IF Iquit=1 THEN 4904 

492@ GOTO 4480 

4904 PRINT "PU SPQ" 

4908 SUBEND 

4912 SUB Stats 

4916 PRINTER IS 70! 


4920 J=@ 
4974 K=Q 
4902. “SEEP 


4952 INPUT “PET FILE 10) ANRSLEYZE?” Files 
4936 ASSIGN @File TO Files 

494Q@ BEEP 

4944 INPUT “LAST RUN No?(@=QUIT)",Nn 
4948 IF Nn=@ THEN S092 
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495. 
4956 
4962 
4364 
4368 
ASi2 
4976 
4980 
4984 
4988 
4992 
4996 
5000 
5004 
5808 
5012 
5016 
50°20 
5024 
5028 
5032 
5036 
5042 
5044 
5048 
5@52 
5QS56 
5062 
5064 
5@68 
SO72 
5@76 


52082 
5084 


Nn=Nn-J 

Sx =Q 

Sy=2 

$2=@ 

Sxs=@ 

Sys=Q 

Szs=Q 

FOR I=1 TO Nn 

J=J+1 

ENTER ®File:Q,T 

H=Q/T 

S~™=Sx+Q 

Sxs=SxstQ*2 

Sy=Sytf 

Sys=SystT*2 

$2=Sz2+H 

Szs=Szst+H*2 

NEXT I 

Qave=S~/Nn 

Tave=Sy/Nn 

Have=5z/Nn 

Sdevqg=SQR(ABS( (Nn#*Sxs-S™*2)/(Nn#*(Nn-1)))) 
Sdevt=SQR(ABS( (Nn#*Sys-Sy*2)/(Nn#*(Nn-l)))) 
Sdevh=SQR(ABS( (Nn#*S28-S272 )/(Nn#*(Nn-l)))) 
Sh=100+*Sdevh/Have 

$Q=100*Sdeva/Qave 

St=100@+Sdevt/Tave 

IF K=1 THEN 5084 


PRINT 

PRINT-USING “11%, “DATA FILE:* (14 ce rres 

PRINT 

PRINT USING “11X,""RUN Htube SdevH Qdp SdevQ Thetab SdevT"“" 
K=1 


PRINT USING *11%,00,6(2%,0.30E,1X%,30.20),2%,00.30,1% 30.20" J, Have jsmecaee 


sods | ave sot 


5288 
5892 
5@96 
5122 
5104 
5188 
ap 
olIT6 
S120 
5124 
5128 
aloe 
5136 
514@ 
5144 
$148 
Sis< 
5156 
5162 
5164 
S168 
Sire 
5175 
518@ 


GOTO 4942 

ASSIGN @Filel TO # 

PRINTER IS 1 

SUBEND 

SUB Coef 

COM /Coly/ AC10,10),C(10),8(4),N,Iornt ,Opo,Ilog 
BEER 

INPUT "GIVE A NAME FOR CROSS-PLOT FILE” ,Cofs 
BEEP 

INPUT “OUTPUT TYPE (@=q vs Dt, l=h vs Dt, 2=h vs gq)" ,Opo 
CREATE BDAT Cof$,& 

ASSIGN @File TO Coffs 


BEER 

INPUT “SELECT (O=LINEAR T=U0G(% ,¥)2",1 hog 
BEEF 

INPUM- "ENTER O1L PERCENT <=14S70P)..86p 
BEEP 


INPUT “ENTER TUBE NUMBER (1, 2, 5, 4, OR S?",Itn 
IF Bop«@ THEN 5176 

Chee eolyGl tn) 

OUTPUT @FileiBop,6/( +) 

GOTO 5144 

ASSIGN @File TO * 

SUBEND 
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5134) 


S188'XXXXXKXXXKXXXXKXKXXXKXKKXKXKXXKX KK KKK KK KKK KK KKK KX KX KXKXKK KK KKK KKK KK 


aes a 
S156 
5220 
5204 
5208 
Se l2 
S2i6 
S22 
5224 
S228 
Gee 
S236 
5240 
$244 
5248 
Syieey 
S226 
5260 
$264 
5268 
Sia i4 
S276 
5280 
5284 
5288 
S252 
S296 
530@ 
5304 
5308 
Sale 
5516 
5320 
5324 
5328 
Boose 
2200 
$340 
5344 
5348 
yeas 
Ses] 
5362 
5364 
5368 
Sere 
b2?6 
5380 
53584 
5388 
So a2 
Bhs |5' 
5400 
5404 
5408 
S412 
5416 
5428 


SUB Plot 
COM /Cply/ AC10,10),C(10),8(4),Nop,Iprnt ,Opec,lIlog 
GIN Xy(17) 
INTEGER [1 
PRINTER IS 1 
BEEP 
Idv=1 
INPUT “LIKE OEFAULT VALUES FOR PLOT (1=Y(OEFAULT),O=N)?7",Idv 
Opo=0 
BEEP 
PRINT USING “4xX,""Select Option: ""“ 
PRINT USING “6X,""0 q versus delta-T(DEFAULT)""" 
PRINT USING "6X,""1 h versus delta-T""" 
PRINT USING "6X,“"2 h versus q*"" 
INPUT Qpo 
BEEP 
INPUT "SELECT UNITS (@=SI¢( DEFAULT), 1=ENGLISH)",Iun 
PRINTER IS 705 
IF Idv<>l THEN 
BEEF 
INPUT “ENTER NUMBER OF CYCLES FOR X-AXIS" ,C~ 
BEEP 
INPUT “ENTER NUMBER OF CYCLES FOR Y=-AxIS* ,Cy 
BEEP 
INPUT “ENTER MIN X-VALUE (MULTIPLE OF 10)" ,Xmin 
BEEF 
INPUT "ENTER MIN Y-VALUE (MULTIPLE OF 10)",Ymin 
a Ss 
IF Opo=@ THEN 
Cy=2 
Cx22 
Xmin=l 
Ymin=1000 
ENO IF 
IF Opo=l THEN 
Cy=2 
Cx=2 
Xmin=] 
Ymin=1Q@ 
ENO IF 
IF OQpo=2 THEN 
Cy=2 
Cx~=2 
Xmin=1000 
Ymin=100 
ENO IF 
ENO IF 
BEEP 
PRINT “IN: SPL. IP 230@,2200,8320,6800: ” 
PRINT "SC 0,100,0,100;TL 2,0: " 
Sf x=100/C.. 
Sfy=100/Cy 
BEEP 
INPUT "WANT TO BY-PASS CAGE (1=Y, O=NO(DEFAULT)",Ibp 
IF Ibp=1 THEN 5908 
PRINT “PU 0,0 PD" 
Nn=9 
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S228. /POR Sa) ies! 


S428 Xatexmine1@*¢( I-1) 

54352 IF T=Cxt+i THEN Nn=1 

5436 FOR J=1 TO Nn 

S44Q IF J=i1 THEN PRINT “TL 2 Q" 
5444 IF J=2 THEN PRINT “TL 1 Q' 
5448 Xa=XateJ 

5452 X=LOT(Xa/Xmin)#Sfe 

S456 PRINT “PA";:X,",@; XTi" 
5460 NEXT) cd 


5464 NEXT I 
S468 PRINT “PA 1@0,@;PU;” 
S472 -PRING PUPA 2,0 Fo- 


5476 Nn=9 

S48@ FOR I=! TO Cy+i 

5484 Yat=Yminel@°¢( I-21) 

S488 IF [=Cy+l THEN Nn=l 

5492 FOR J=1 TO Nn 

5496 IF J=l1 THEN PRINT “TL 2 Q” 
5500 IF J=2 THEN PRINT "TL 1 Q" 
5504 Ya=Yat#J 

5508 Y=LGT(Ya/Ymin)#Sfy 

S)=)) 74 PRINT "PA 0,°;:Y,"YT* 

S516 NEXT J 


S552@ NEXT I 
5524 PRINT “PA 0,100 die @ 2° 


S528 Nn=9 

B532 FOR vI=) 10) Cx+1 

S536 Xat=Xmin#l@*cI-1) 

5540 IF I=Cx+i1 THEN Nn=} 

5544 FOR J=1 TO Nn 

5548 IF J=1 THEN PRINT “TL @ 2" 
6552 IF J>1 THEN PRINT “TL @ 1" 
5556 Xa=XateJ 

5560 X=LGOT(Xa/Xmin)#Sfx 

5564 PRINT “FA-ixk; ,10G; x1 - 
S568 NEXT J 


S572 NEXTol 
S5/G PRINT “PA 100,100 PU PA 100,060 FD™ 


558@ Nn=3 

5584 FOR [I=1 10 Cy+!1 

5588 Yat=Yminel@*(I-1) 

SSig2 IF T=Cy+i THEN Nn=1l 

5596 FOR J=1 TO Nn 

5600 IF J=l1 THEN PRINT "TL @ 2" 
5604 IF J>1 THEN PRINT "TL @ 1" 
56@8 Ya=YateJ 

S612 Y=L6T(Ya/Ymin)*Sfy 

5616 PRINT “PO PA 100,",Y, "YT" 
5622 NEXT J 


S624 NEXT I 

S628 PRINT “PA 100,100. Fu: 
S632 PRINT “PA @,-2 SR !1.5,2" 
S636 J1=LOT( mind 

5640 FOR IT=1 TO C.+1 


5644 ~ XKa=XminelO°CI-1) 

5648 X=LGT(XKa/Xmind#Sf« 

5652 PRIN eA eee oe 

S656 IF Ii1e#@ THEN PRINT "CP -2,-2;LB1@;:PR -2,2;LB":11;"" 
5660 IF I1-@ THEN PRINT “CP -2,-2:;LB10;PR 0,2:LB°:32.°"" 
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S664 Ti=ii+l 

S668 NEXT I 

S672 PRINT “PU PA 2,0" 
S676 [i=LGT(Ymin) 
S680 Y10=10 

S684 FOR I=1 TO Cytl 


5688 Ya=Ymin*l@*(I-1) 

5692 Y=LGT(Ya/Y¥min)#Sfy 

S696 PRINT "PA 0,°;:Y,°" 

S702 PRINT “CP .=4..—.25:LB10;PR -2,2:LB°:311,"" 
5724 Ti=lIitl 

S708 NEXT I 

Sel2 “BEEP 

S716 Idle=l 


S720 INPUT “WANT USE DEFAULT LABELS (1=Y(DEFAULT),@=N)7",Id1 
$724 IF Idil<>t THEN 


5728 SEE 
S732 INPUT "ENTER X-LABEL",XlabelS 
$736 BEEP 
5742 INPUT "ENTER Y-LABEL",Ylabel$ 


S744 ENO IF 
S748 IF Qpo<2 THEN 


S7S2 PRINT "SR 1,2;PU PA 40,-14; " 

S756 PRINY “EBCT. ER =126.5 PO PR 1.250 PU;PR .5.-4; LBuwo;PR .5,1;" 

5762 PRINT LB=14PR 25, -Il;LGsat>PR <5, 13° 

S764 IF Iun=@ THEN 

S768 PRINT “EB} / (K}° 

6772 ELSE 

5776 PRINT "LB) / (F)" 

5780 END IF 

5784 END IF 

S788 IF Opo=2 THEN 

S792 IF Tun=@ THEN 

S796 PRINT "SR 1.5,2;PU PA 40,-14;:LBq / (W/m;SR 1,1.5;PR @.5,1;LB2;SR 1 

Poacite O.5,-1;3;LB)~ 

5802 ELSE 

5804 PRINT “SR 1.5,52;PU.Fe 34, -14;LBa / 4Btu/hrs PR .5..5¢LB.;PR .5,=.5: 

5808 PRINT “UBPFtrPR..Sel. oR 121.5:ELB2:SR 1.5, 23Pk.S5,-1;3L8):” 

5812 END IF 

S816 ENO IF 

5820 IF Opo=@ THEN 

$824 IF Tun2@ THEN 

S828 PRINI “SR 1.5,2:3FU PA =12 40,01 0, 1;LBa 7 (W/m:cPR -1,8-.5;: SR 11.5;L 

B2;SR 1.5,2:PR 1,.5;LB)" 

Bese ELSE 

$836 PRINT "SR 1.5,2;PU PA -12,52:01 0,1:LBq /¢ (Btu/hr;PR -.S,.5:LB.:PR 
2 eae 

5842 PRINT “L6ft;SR 1,1.S:;PR -1,.S;LB2;PR 1,.5,SR 1.5,2;:LB)" 

S44 ENO IF 


S848 END IF 
Sage IF Opo*0 THEN 


S856 IF Iun=@ THEN 

5862 PRIN Ase 125 2¢FU PA 12 oe blo. tLe hay rum PR = 1 6 Sask 1.1.5. 2S 

225k 1.5,2+;PR-.5,457 " 

S864 . PRINT “LB.: PR .S,@:LBK)" 

S868 alse 

5872 PRINT "SR 1.5,c2:iPU PA -12,28:01 @,1:LBh / (Btu/hriPR -.S,.5;LB.:PR 
Sans 

5876 PRINT MOUS Cte i. o.5r 1.15: UR SSR See sPR 5, . Shee G PR 9.5.4. 5% 
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LBF >” 


Sa20 END IF 

S884 END IF 

S888 IF Idl=@ THEN 

S892 PRINT “SF 1.S,2:PU PA SO,-16 CP”: —-LEN(Xlabel$)/2; "0;LB": xX labels ae 
S896 PRINT "PA -14,50 CP 0,";-LEN( Ylabel$)/2+*5/6; "DI 0,1:LB": Ylabel$;i ~~ 
5980 PRINT “CP @,@ DI" 

S904 END IF 

5908 Ipn=0 

S912 Repeat: ! 

S916 Xll=!.E+6 

S920 Xul=-1.6+6 

5924 Iecn=0 

S928 BEEP 

S35. Sela 

S936 INPUT "WANT TO PLOT DATA FROM A FILE (1=Y(DEFAULT),@=N)?" ,Ok 
S940 IF Ok=1 THEN 

5944 BEEP 

S948 INPUT “ENTER THE NAME OF THE DATA FILE” ,OfileS(Q) 
S9S2 ASSIGN @File TO OfileS(Q) 

S956 BEEP 

S96@ Npairs=20 

S964 INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED( DEFAULT=20)",Npairs 
S968 BEEP 

S972) [Itn=Itnel 

S976 INPUT "ENTER TUBE NUMBER (1, 2, 3, 4, OR S)",Itn 
5980 BEEP: 

5984 PRINTER IS 1 

5988 INPUT “WANT DEFAULT SYMBOLS? (YES=@ (DEFAULT),NO=1)" ,Symb 
5992 Sym=Itnt2 

S996 IF Symb=Q@ THEN 

6800 GOTO 6036 

6004 END IF 

60@8 PRINT USING "4X,""Select a symbol:""" 

6012 PRINT USING "6X,°"1 Ster 2 Plus sign""" 
6016 PRINT USING "6X,""3 Circle 4 Square""" 

6020 PRINT USING "6X,""S  Rombus""" 

6024 PRINT USING "“6X,""G Right-side-up triangle""" 
6028 PRINT USING "6X,""7 Up-side-down triangle"”"" 
6@32 INPUT Sym 

6036 PRINTER IS 70S 

6@4@ PRINT “PU OI" 

6044 IF Sym=! THEN PRINT "SM?" 

6048 IF Sym=2 THEN PRINT “SM+t" 

6@S2 IF Sym=3 THEN PRINT “SMo" 

6056 FOR I=! TO Npairs 

6@6@ ENTER @Fileixy(*) 

6064 IF Spo=@ THEN 

6068 Ya=Xy(Itn-1) 

6072 Xa=Xy(lil¢+Itn? 

6@77 END IF 

6@8@ IF Opo=!l THEN 

6@84 Ya=xy(StItn) 

6058 Ka=xy(lit+Itn>? 

6@92 END IF 

6096 IF’ Opo=2 THEN 

6100 Ya=Xy(S+Itn) 

6104 Xa=Xy(Itn-1) 

6108 END IF 

6Ette IF XasXll THEN X11=Xa 
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Bllo IF Xa Kul THEN Xul=Xa 
612@ IF TIun=1 THEN 


6124 IF Opo<2 THEN Xa=Xael.8 

6128 IF Opo>@ THEN Ya=Ya*.1761 
6132 IF Opo=@ THEN Ya=Ya*.317 
6136 IF Opo=2 THEN Xa=Xae.317 


6140 END IF 

6144 X=LGT(Xa/Xmin)+#Sfx 

6148 Y=LGT(Ya/Ymin)+#Sfy 

6152 kK j)=0 

61S6 CALL Symb(X,Y,Sym,Icl,Kj) 

6160 GOTO 6212 

6164 IF Sym>3 THEN PRINT “SM° 

6168 IF Sym<4 THEN PRINT "SR 1.4,2.4" 

6172 IF Icl=@ THEN 

6176 PRINT “PA",X,Y,"" 

6180 ELSE 

6184 PRINT “PA’.%,Y, PD" 

6188 ENO IF 

6192 IF Sym>3 THEN PRINT "SR 1.2,1.6" 

6196 IF Sym=4 THEN PRINT "UC2,4,99,0,-8,-4,90,9,8,4,@: " 
6200 IF SymeS THEN PRINT *°UC3,0,99,-3,-6,-3,6,3,5,3,-5: ° 
B2045 IF Sym=6 THEN PRINI "UCO,5.3,53,5,-9,;-6,0,5.35: © 
6208 IF Sym=7 THEN PRINT “UCO,-5.3,99,-3,8,6,0,-3,-8: “ 
6212 NEXT I 

6216 PRINT "PU" 

6220 BEEP 

6224 Ilab=l 

6228 INPUT "WANT TO LABEL? (1=Y(DEFAULT),Q@=N)",Ilab 
6232 IF Ilab=1 THEN 


6236 PRINT "“SPQ;SP2” 

6240 BEEF 

6244 IF Klab=@ THEN 

6248 Xlab=65 

6252 Ylab=85 

6256 INPUT “ENTER INITIAL X,Y LOCATIONS" ,Xlab,Ylab 
6262 Xtt=Xlab-S 

6264 Ytt=Ylab+8 

6268 PRINT "SR 1,1.5" 

6272 PRINT “SM:PA" ,Xtt,Ytt, "LB Tube y File" 
6276 Ytt=Ytt-3 

6282 PRINT “PA” ,Xtt,Ytt, “LB No O1l Name" 
6284 IF Sym=1 THEN PRINT “SMe” 

6288 IF Sym=2 THEN PRINT "SM+" 

Bede IF Sym=3 THEN PRINT "SMo" 

6296 Klab=1 

6300 END IF 

6304 kKy=1] 

6308 CALL Symb(Xlab,Ylab,Sym,Icl,k3) 

b3l2Z PRINT “SR 1,1.5:;SM”" 

6316 IF Sym: 4 THEN PRINT “PR 2,0" 

5329 PRINT “PR 2,0;L6 eling 7 

6324 BEEP 

6328 INPUT “ENTER BOP‘ Q=DEFAULT)" ,Bop 

b332 IF Bop<1@ THEN PRINT "PR 3,0:LB"i Bop: “" 

6336 wile (Boo7S THEN PRINT. "PR 125.0706" *S6p: ~~ 

6340 PRINT “PR 2,0:LB";:OfilesS(Q@)i"" 

6344 PRINT "SP@:SP1;SR 1.5,2” 

6348 Ylab=Ylab-S 


6352 END IF 
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a a 
ASSIGN ®6File TO » 
XLL=X11/1.2 
Xul=Xul#l.2 
GOTO 8042 
END IF 
PRINT “PU SM" 
BEEP 
Go_on=1 
INPUT “WANT TO PLOT A POLYNOMIAL (1=Y(OEFAULT ),Q0=N)?",Go_on 
IF Go_on=1 THEN 
BEEP 
PRINTER IS 1 
INPUT "WANT DEFAULT LINE TYPE? (YES=@ (DEFAULT),NO#=1)",Ln 
Ipn=Itn 
IF Ln=@ THEN 
GOTO 6448 
END IF 
PRINT USING "4X,""Select line type:""" 
PRINT USING "6X,""@ Solid line. * 
PRINT USING "6X,""1 Dashed"“" 
PRINT USING "“6X,""2,,,5 Longer line - dash""" 
INPUT Ipn 
PRINTER IS 7@S 
BEEP 
Ilog=1 
INPUT "SELECT (Q@=LIN,1=2LOG(DEFAULT))" ,Ilog 
Iprnt=l 
CALL Poly(Ofile$(*),Npairs,Itn) 
FOR X~=@ TO Cx STEP Cx/2Q@ 
Xa=XminelO*xX-« 
IF Xa<Xll OR Xa>Xul THEN 6572 
Icn=Icntl 
Pu=@ 
IF Ipn=1 THEN Idf=Icn MOD 2 
IF Ipn=2 THEN Idf=Icn MOD 4 
IF Ipn=3 THEN Idf=Icn MOD 8 
IF Ipn=4 THEN Idf=Icn MOD 16 
IF Ipn=S THEN Idf=Icn MOD 28 
IF Idf=1 THEN Pu=l 
Ya=FNPoly( Xa) 
IF YacYmin THEN 6572 
IF Iun=1 THEN 
IF Opo<2 THEN Xa=Xael.8 
IF Opo>@ THEN Ya=Ya*.1761 
IF Opo=@ THEN Ya=Ya*.317 
IF Opo=2 THEN Xa=Xa*.317 
END IF 
Y=L6T( Ya/Ymin)#Sfy 
K=LoTixa/Xmind#Sfa« 
IF ¥:@ THEN Y=Q 
IF Y:10@ THEN GOTO 65972 
IF Pu=@ THEN PRINT “PA",xX,Y, "PD" 
IF Pu=l THEN PRINT “PA",X,Y, "PU" 
NEXT X. 
PRINT "PU" 
ENO IF 
BEEP 
INPUT "WANT TO QUIT (1=Y,@=N(DEFAULT>)",Iat 
IF Iqt=1 THEN 6602 


LiZ 


Boob couUle-agic 

66@@ PRINT "FU PA 8,0 SPQ" 

66@4 SUBEND 

66Q8 |! 
BEL2IXXXXXXKXXXXXXXXXXXAXXXKKAXAKXXKXXXXKXXKKKXXKXXAKXKKXKKKXKKK KK KKK XXKXKKKKK KK 
6616! 

6620 SUB Symb(X,Y,Sym,Icl,Kj) 

6624 IF Sym>3 THEN PRINT “SM” 

6628 IF Sym<«4 THEN PRINT “SR 1.4,2.4" 
6632 Yad=0 

6636 IF Kj=1 THEN Yad=.8 

664@ IF Icl=@ THEN 


6644 PRINS cPhra kA, YryYad, -— 
be45 ELSE 
ies BRIN “PA ox ,\tvYad, PD 


6656 END IF 

6660 IF Sym-3 THEN PRINT “SR 1.2,1.6” 

6664 IF Sym=4 THEN PRINT "UC2,4,99,0,-8,-4,@,0,8,4,Q0;" 
6668 IF Sym=S THEN PRINT "“UC3,0,99,-3,-6,-3,6,3,6,3,-6." 
6672 IF Sym=6 THEN PRINT “UCO,5.3,99,3,-8,-6,0,3,8:;" 
Be7G ste Syme? THEN PRINT “UCO,-5.3,99,-3.8,6,0,-3.-8: ’ 
6680 IF Kjy=1 THEN PRINT “SM:PR @,~.8" 

6684 SUBEND 

6688! 

6692 '!XXXXXXXXXXXXXXXXXXXXXXXXXKXKXXKXKXKXXKKXXKKXKXK XX KK KK KXKXKKXKXKXKKKXKKXKKKKKK KKKK 
6696! 

67@0@ SUB Fixup 

B7g4! FILE: FIXUP 


6708! 
6712 DIM Emf(34),Amp(11),Volt(4),Ldtc(4) 
p¢iG. BEEP 


G7¢2O INPUT “OLD FILE TO FIXUP” O2files 

6724 ASSIGN @File2 TO D2files 

6728 OlfilesS="TEST" 

6732 CREATE BDAT Difile$,6@ 

6736 ASSIGN ®Filel TO Difiles 

674@ ENTER ®File?:Nrun,Date$,Ldtc(+),Itt,Bop,Nht ,Natp,Nrt ,Corr 
6744 OUTPUT ®Filel:iNrun,Date$,Ldtc(*),Itt,Bop,Nht ,Natp,Nrt,Corr 
6748 FOR I=1 TO Nrun 


67582 ENTER ®File2; Told$,Emf(+*),Volt(#),Amp(#) 
6756 IF I=] THEN 6764 
676@ OUTPUT @Filel: Bop, ToldS,Emf(#),Volt(*),Amp( +) 


6764 NEXT I 

6768 ASSIGN @File2 TO * 

6772 ASSIGN OFilel TO » 

Bava wRENAME “TEST” TO D2 files 
Bie@ *csEEP [000,. 
Bvs4 BEEP 4000,. 
6788 BEEF 400@,. 
e7ge.  SuUBEND 
6796! 

BSA! KKM KM KM MK KKM KKK MMM KK KX KKM KKM KKK KK MM MK KRM XK KK OK KK KY 

6204) 

68@8 SUE Move 

B812! FILE NAME. MOVE 

6816! : 

6820 DIM A’ 66,,B8(66),C(66),0(66),E(66),F:56),6(66),H(65),3(66),k(66?,L166),M(66 
) 

pact PIM N( G6) Emfi 54), Volt(2),AmeC ll .bdte(4) 

soce BEEP 


eee a Pl 8 | 
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bd5 2. INU a ORDSF LER TOsMOUES De i1 es 

6856 ASSIGN @Filee TO O2_files 

684@ ENTER @File2iNrun,Told$,Ldtc(*),Itt,Bop,Nht,Natp,Nrt,Corr 
6844 FOR [=1 TO Nrun 


6848 ENTER ®File2:iTold$ 

6852 ENTER @FileziAC lI), ,6(1),C¢(1),0(1),6¢(1),F¢(1),601 ) Bal), 1) Kt 1) eae 
I),NC(I) 

6856 ENTER @File2:Emf(*#),Volt(#),Amo(*#) 


6860 NEXT I 

6864 ASSIGN ®6File2 TO * 

6868 BEEP 

6872 INPUT “SHIFT DISK AND HIT CONTINUE”, Ok 

6876 BEEP 

688Q INPUT "INPUT BDAT SIZE",Size 

6884 CREATE BDAT D2_file$,Size 

6888 ASSIGN @Filel TO D2_files 

6892 OUTPUT @FileliNrun,Date$,Ldtc(*),Itt,Bop,Nht,Nato,Nrt,Corr 
6896 FOR I=1 TO Nrun 


6922 OUTPUT @FileliTolds 

6904 OUTPUT GFilel: ACT) BCI) CCl), O(1) ECl) Fi 1) 6( 1 Att), 361) Re ee 
CT) NK) 

6988 OUTPUT @FileliEmf(#),Volt(#*),Amp(*) 


6912 NEXT I 

6916 ASSIGN @Filel TO * 

6920! RENAME “TEST 1G DO2_f mies 

6924 BEEP 2000,.2 

6928 BEEP 4000,.2 

B9S2556EEP 4000.2 

6936 PRINT “DATA FILE MOVED” 

694Q@ SUBEND 

6944! 

6948 '!XXXXXXXXXXXXKXKXXXKXXXKXKXKXKXKXXKXKXKXKXKXKXKXKXKXXKKXKKXKKKKXKKKKKKKKK KKK KK KKK KK KX 
6952! 

SIS56085U6 Pom 

69650 8EEP 

6964 INPUT "ENTER FILE NAME TO BE DELETED” ,Files 
6968 PURGE Files 

6972 GOTO 6960 

6976 SUBEND 

698Q! 

6984! XXXXXXXXXXXXXXXXXXKXAXAXAXXKXKXAXKXXXKKXAXXKXAKXAXKAKAKXKXKX KKK XA AKKXKXKAKXXXKXAKXKXKX KKK KX XKAK 
6988! 

6gg25° SUB Come 

6996! FILE NAME: COMB 


7QQQ ! 
7004” DIM Emf(54),Volt(Z),Ampt tl? batet4) 
7@08 GEEP 


(Olle INPUT “OLOSFILE TO FIXUP {Dc files 

7016 ASSIGN @Fiie2 TO De_files 

7020 Ol_fileS="TEST" 

7024 CREATE BDAT DI_files$,3@ 

7028 ASSIGN @File! TO DI_files 

7032 ENTER GFilel:Nrun,Oate$,Ldtc(*),I[tt,Boo,Nht,Nato,Nrt,Corr 
7036 IF K=@ THEN OUTPUT @FileliNrun,Date$,idtc(*#),Itt,Bop,Nht ,Nato ,Nrt, Corr 
7040 FOR I=! 1S Nrun 

7244 ENTER @File2:Bop,ToldS,Emf(+«) Volt(+),Amp(e) 

7048 OUTPUT @Filel:Bop,ToldS,Emf(*),Vclt(*#),Amp(#) 

70@S2 NEXT I 

7@S6 ASSIGN @Filec TO * 

706@!' RENAME "TEST" TO Ov_files 
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7@E4 
7068 
7072 
7876 
708 
7084 
7088 
70392 
7096 
7100 
7104 
7108 
tide 
7116! 


BEEP 4000,.2: 

BEEP 

Oka=1 

INPUT "WANT TO ADD ANOTHER FILE (1=#Y,Q@=N( default ))?7",Oka 
IF Oka=1 THEN 

K=1 

SEES 

INPUT "GIVE NEW FILE NAME" ,Nfiles 
ASSIGN @File2 TO Nfiles 

S0TO 7032 

END IF 

ASSIGN ®Filec TO * 

SUBEND 


TILQEXXXXXKXKXKXXXXXXXXKXAXXAXXXKXKXAXAXKXK AX KKK KAKAKKKXXAXKAKXAKXKXKXKXAAXAXXKXKAKXXKXXXAKAKAKK 


7124! 
T1258 
(er 
7136 
7148 
7144 
7148 
7152 
7 256 
7160 
7164 
7168 
mate 
7176 
718 
7184 
7188 
Tide 


SUB Readplot 
DIM Qdp(5),Htube(5), Thetab(S) 
PRINTER IS 701 
INPUT "ENTER FILE NAME" ,File$ 
INPUT “ENTER THE NUMBER OF DATA PAIRS" ,Neun 
ASSIGN @Filel TO Files 
FOR I=] TO Nrun 
ENTER @Filel:iQdp(«),Htube(#),Thetab(«) 
PRINT Qdp(*) 
PRINT 
PRINT Htube(#) 
PRINT 
PRINT Thetab( #) 
PRINT 
PRINT 
NEXT I 
SUBEND 
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